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A  BRIEF  SIMMARY  OF  THE  RESEARCH  ACHIEVEMENTS  DURING 

CY  1977-1981 

At  the  beginning  of  1977,  the  possible  relationship  between  the  transition¬ 
al  spot  and  the  large  coherent  structure  in  a  fully  developed,  turbulent 
boundary-layer  was  recognized.  A  portion  of  the  spot  seemed  to  retain  its 
identity  without  observable  deterioration  for  at  least  100  boundary-layer 
thicknesses  (Zilberman,  Wygnanski ,  Kaplan,  1977).  The  experiment  was  repe¬ 
ated  at  USC  for  the  purpose  of  verification  and  a  possible  sharper  defini¬ 
tion  of  the  scales  in  the  spanwise  direction.  The  second  experiment  incor¬ 
porated  two  detector  probes  used  for  selection  of  eddies  which  were  sym¬ 
metrically  aligned  with  an  array  of  hot  wires.  The  latter  recorded  the  en¬ 
tire  velocity  profile  at  any  Instant  rather  than  the  velocity  at  a  single 
point,  which  was  recorded  by  Zilberman,  et.  al.  1977.  The  finds  were 
quite  similar  (Harltonldis,  Kaplan,  Wygnanski,  1978)  to  the  ones  esta¬ 
blished  earlier.  One  could  evaluate,  however,  the  effects  of  signal  condi¬ 
tioning  on  the  flow  field  educed,  because  instantaneous  velocity  data  span¬ 
ning  the  entire  boundary  layer  was  available. 

The  difference  between  the  scale  of  the  transitional  spot  embedded  in  a 
laminar  boundary  layer  and  the  scale  of  the  structure  educed  in  a  fully 
turbulent  environment  led  to  many  pressing  questions,  some  of  which  are 
listed  below: 

(1)  Is  the  spot  a  single  large  coherent  eddy  on  whioh  small  scale 
turbulence  is  superimposed  or  is  it  an  assembly  of  eddies,  both 

large  and  small? 
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(2)  How  does  the  spot  depend  on  the  character  of  the  surrounding 
boundary  layer? 

(3)  How  does  the  spot  grow?  Does  it  entrain  external  fluid  or  does 
it  destabilize  the  surrounding  boundary  layer  -  or  both? 

(4)  Is  there  any  way  by  which  the  rate  of  growth  of  the  spot  can  be 
controlled? 

(5)  How  do  spots  behave  when  they  interact  with  one  another,  either 
longitudinally  or  laterally? 

(6)  What  are  the  effects  of  Reynolds  number  and  pressure  gradient  on 
the  rate  of  growth  of  the  spot? 

Some  of  the  recent  investigations  were  aimed  at  answering  these  questions. 

An  oblique  packet  was  found  to  trail  a  transitional  spot  (Wygnanski,  Hari- 
tonidis,  Kaplan  1979)  at  its  wingtlps.  The  maximum  spreading  angle  of  the 
packet  is  11°,  which  is  identical  to  the  lateral  spread  of  the  spot,  sug¬ 
gesting  a  possible  relationship  between  the  two.  The  packet  was  related  to 
the  Tollmlen-Schllchtlng  wave  instability  and  the  waves  detected  underwent 
the  strongest  amplification,  according  to  the  linear  stability  of  the  Hla- 
sius  velocity  profile.  At  higher  Reynolds  number  (Re),  these  waves 
broke-down  and  generated  new  spots,  which  trailed  the  original  spot  near 
its  wing-tips. 
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The  presence  of  the  wave  packet  and  its  breakdown  suggests  that  the  spot  is 
not  a  single  coherent  eddy,  but  rather  an  orderly  formation  of  large  ed¬ 
dies.  Furthermore,  the  lateral  spread  of  the  spot  may  have  been  caused  by 
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the  destabilization  of  the  laminar  boundary  layer.  These  finds  provided 
the  first  link  between  stability  theory  and  the  actual  spreading  of  turbu¬ 
lence.  We  expected  the  Reynolds  number  to  play  an  important  role  in  the 
spreading  of  turbulence  in  laminar  boundary  layer,  because  the  instability 
of  the  Blasius  boundary  layer  depends  on  viscosity.  There  was,  however,  no 
data  available  which  showed  any  dependence  of  the  spot  on  Re.  A  systematic 
study  (Wygnanski,  Haritonitis,  Zilberman  19S0)  showed  that  the  rate  of 
elongation  of  the  spot  in  the  streamwise  direction  is  dependent  on  the  Rey¬ 
nolds  number  existing  at  the  origin  of  the  spot.  The  differences  in  the 
rate  of  elongation  of  the  spot  are  limited  to  the  celerity  of  the  trailing 
turbulent-non-turbulent  interface  of  the  spot  which  decelerates  with  incre¬ 
asing  Re.  It  appears,  therefore,  that  the  trailing  edge  of  the  spot  is 
probably  most  active  in  destabilizing  the  boundary  layer.  The  spreading 
rate  of  the  spot  in  the  spanwise  direction  is  independent  of  Re  while  the 
height  of  the  spot  increases  with  downstream  distance  in  the  same  manner  as 
the  turbulent  boundary  layer  does,  and  is,  thus,  only  weakly  dependent  on 
Re.  The  different  functional  dependence  of  the  growth  of  the  spot  in  the 
three  principal  directions  on  Re  implies  that  the  spot  does  not  obey  a  un¬ 
iversal  conical  similarity  proposed  by  Cantwell,  Coles  and  Dimotakis 
(1978);  it  also  follows  that  there  is  no  single  mechanism  responsible  for 
the  spreading  rate  of  the  spot.  Thus,  while  the  aft  of  the  spot  is  associ¬ 
ated  with  the  breakdown  of  Tollraien  Schlichting  waves,  the  center  and  the 
front  of  the  spot  may  simply  engulf  and  entrain  non-turbulent  fluid 
Detailed  calculations  of  entrainment  into  an  ensemble  averaged  boundary  of 
the  spot  on  its  plane  of  synmetry  were  made  by  Zilberman  (1982),  who  used  a 
rake  of  x-wire  anemometers  in  his  calculations. 
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This  investigation  revealed  that  entrainment  calculations  based  on  ensemble 
averaged  data  in  the  x-y  plane  are  Inadequate  because  significant  entrain¬ 
ment  must  occur  in  the  x-z  plane.  Velocity  contours  reconstructed  from  in¬ 
dividual  realizations  also  suggested  that  the  spot  contains  numerous  eddies 
which  engulf  non-turbulent  fluid. 


Favorable  pressure  gradient  decreased  significantly  the  rate  of  elongation 

of  the  spot  which  dropped  from  =  0.45  to  0.29  (L  being  the  length  of 

dx 

the  spot)  at  a  comparable  Reynolds  number.  The  spanwise  spreading  angle  of 
the  spot  was  reduced  from  11°  to  5°.  These  significant  differences  in  the 
contamination  of  the  laminar  boundary  layer  by  turbulence  occurred  in  spite 
of  the  fact  that  the  effect  of  favorable  pressure  gradient  on  the  laminar 
boundary  layer  profile  was  rather  small.  (The  measured  profile  correspond¬ 
ed  to  a  Falkner-Skain  profile  with  a  parameter  =  0.12).  However,  very 
slight  favorable  pressure  gradient  is  required  in  order  to  appreciably  im¬ 
prove  the  stability  of  the  boundary  layer.  The  slower  destabilization  of 
the  laminar  boundary  layer  enabled  us  to  detect  the  existence  spanwise  in- 
hcmogeneties.  which  appeared  even  in  a  simple  ensemble  averaged  data.  A 


preselection  of  the  events  before  averaging  revealed  clearly  the  existence 
of  longitudinal  eddies  in  the  spot.  It  was  the  first  proof  that  the  spot 
consists  of  an  array  of  large  hairpin  eddies  moving  downstream  in  a 
formation.  Measurements  at  various  downstream  distances  revealed  that 
the  r.unber  of  the  hairpin  eddies  across  the  span  of  the  spot  increased  with 
(Wygr.anski  1980)  in  the  number  of  eddies  which  is  probably  caused  by  con¬ 


tinuous  breakdown  of  wave-crests. 
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Measurements  aimed  at  identifying  the  large  coherent  structures  inside  the 
spot,  in  the  absence  of  pressure  gradient,  were  undertaken.  The  Reynolds 
number  in  this  investigation  was  low  in  order  to  detect  the  large  eddies 
more  readily  and  develop  a  rational  scheme  for  the  analysis  of  the  data 
educing  a  most  probably  ensemble  averaged  flow  field,  which  would  be  favor¬ 
ably  compared  with  an  actual  single  realization  (Wygnanski  1981).  This 
scheme  revealed  that  the  number  of  eddies  in  a  spot  increases  linearly  with 
increasing  downstream  distance  in  both  streamwise  and  spanwise  directions. 
The  spot,  thus,  contains  a  large  number  of  coherent  structures  resembling  a 
turbulent  boundary  layer.  What  distinguishes  the  spot  from  a  fully  devel¬ 
oped  turbulent  boundary  layer  is  the  preferred  location  of  the  large  eddies 
with  respect  to  the  spot  boundaries.  Thus,  if  we  are  interested  in  clari¬ 
fying  the  dynamics  of  these  large  eddies  the  transitional  spot  offers  a  un¬ 
ique  possibility. 

The  study  of  transition  in  a  pipe,  which  is  related  to  an  investigation  of 
large  coherent  eddies  existing  in  internal  flows,  proceeded  Intermittently 
over  these  years.  It  was  established  that  the  "slug"  which  is  a  very  large 
turbulent  structure  whose  interior  is  statistically  identical  to  a  fully 
developed  turbulent  pipe  flow,  is  generated  by  an  array  of  "puffs"  merging 
with  one  another.  The  latter  are  orderly  structures  of  clearly  identifi¬ 
able  length  scale  and  character.  The  puff  consists  of  a  small  number  of 
predominantly  toroidal  eddies  which  represent  a  characteristic  transitional 
mode  in  the  pipe.  Periodic  oscillations  were  introduced  into  the  flow  in 
order  to  assess  the  relaxation  times  of  turbulence  and  possible  effect  of 
viscosity  (l.e.,  on  Re)  on  the  large  turbulent  structures. 
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If  I  had  to  sunnarize  our  research  achievements  in  one  sentence,  I  would 
state  that  a  link  has  been  established  between  transition,  on  one  hand,  and 
the  large  coherent  structures  governing  turbulent  shear  flows  at  large  Re 
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abstract 


The  final  steps  of  transition  and  the  generation  of  a  turbulent  boun¬ 
dary  layer  was  simulated  by  generating  a  succession  of  spots  and  allowing 
them  to  interact  under  controlled  environment.  Spots  were  also  generated 
simultaneously  at  two  points  separated  in  the  spanwise  direction  and  their 

mode  of  interaction  was  observed.  Measurements  were  made  on  a  flat  plate 

5 

at  a  representative  Reynolds  number  of  7x10  and  using  two  types  of  rakes 
of  hot  wires:  (i)  wires  normal  to  the  flow  direction  and  parallel  to  the 
surface,  measuring  instantaneously  the  streamwise  component  of  velocity 
across  the  entire  boundary  layer,  (ii)  A  rake  of  wires  which  are  inclined 
to  the  surface  and  to  the  mean  flow,  measuring  instantaneously  the  stream- 
wise  and  the  normal  components  or  velocity  across  the  flow.  The  data  was 
processed  digitally  and  the  information  gathered  was  used  to  assess  the 
statistical  methods  commonly  employed  in  gathering  data  in  a  boundary  layer 
undergoing  transition.  It  is  established  that  ensemble  averaged  velocities 
do  not  represent  the  details  of  the  flow  adequately  and  ensemble  averaged 
turbulent  intensities  may  be  physically  deceiving. 

The  interaction  of  spots  generated  in  succession  from  a  single 
point-source  provides  a  reason* Me  simulation  of  the  turbulent  boundary 
layer.  This  mode  of  interaction  inhibits  the  growth  of  individual  spots  in 
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the  streamwise  and  in  the  spanwise  direction.  As  a  result  of  this  interac¬ 
tion  the  largest  excursions  of  the  uv  product,  representing  the  instantane¬ 
ous  shear  stress,  occur  near  the  leading  interface  of  each  spot.  The  en¬ 
trainment  of  non  turbulent  flow  in  the  X-Y  plane  can  only  account  for  20% 
of  the  total  entrainment  necessary  to  sustain  the  growth  of  the  turbulent 
boundary  layer  which  leads  one  to  suppose  that  the  entrainment  in  the  Y-Z 
plane  resulting  from  vortices  oriented  in  the  streamwise  direction  is  domi¬ 
nant.  Zone  averaged  turbulent  intensities,  Reynolds  stress  and  time-mean 
velocity  profiles  measured  in  a  succession  of  interacting  spots  and  in  a 
fully  developed  turbulent  boundary  layer  are  in  good  agreement  with  one 
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miRQOUCTIQH 

.  There  is  an  abundance  of  evidence  that  the  fluid  notion  within  the 
turbulent  boundary  layer  is  not  as  randan  as  it  was  once  believed,  but  con¬ 
tains  large  coherent  eddies  whioh  nay  govern  sone  important  features  of 
this  flow.  Although  some  characteristic  behaviour  of  these  eddies  was  ob¬ 
served  and  charted,  the  information  available  to  date  is  still  mostly  des¬ 
criptive  in  nature.  There  are  considerable  difficulties  in  defining  a  "de¬ 
tector"  which  is  capable  of  identifying  large  eddy  structures  unambiguously 
ind  impartially.  Some  difficulties  are  attributed  to  the  fact  that  we  are 
concerned  with  repeatable  events  occurring  quite  randomly  in  space  and 
time.  Furthermore,  these  structures  are  embedded  in  an  environment  con¬ 
taining  a  wide  spectrum  of  finer  scales.  Thus  we  cannot  define  a  signature 
of  a  large  eddy  without  a  priori  knowledge  of  its  shape  and  location  rela¬ 
tive  to  the  observer  and  we  can  not  map  such  an  eddy  without  having  a  prop¬ 
er  criterion  for  pattern  recognition. 

Another  difficulty  in  charting  the  large  eddy  stems  from  the  fact  that 
at  any  instant  of  time  it  occupies  a  volume  of  fluid,  while  most  quantita¬ 
tive  measuring  devices  provide  information  at  a  point.  For  these  reasons 
major  contributions  towards  the  recognition  of  coherent  structure  were  made 
by  visual  methods.  Blaokw*»lder  and  Kaplan  (1972)  are  credited  to  be  the 
first  ones  to  use  an  array  of  sensors  spanning  the  boundary  layer  and  pro¬ 
viding  instantaneous  velocity  information  across  the  flow,  but  even  this 
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information  is  insufficient  to  describe  a  three-dimensional  structure  quan¬ 
titatively. 

A  possible  solution  of  the  dilemma  oould  be  provided  by  perturbing 
laminar  flow  from  a  point-source  in  a  way  which  evokes  the  generation  of 
turbulent  spots.  In  this  way  one  obtains  both  a  time  reference  and  a  mean 
trajectory  along  which  the  spot  travels.  Moreover,  because  the  events  are 
repeatable,  the  smaller  scales  which  are  randomly  superimposed  on  them  van¬ 
ish. 

The  difficulty  with  this  approach  in  comparison  with  the  conventional 
"chasing"  of  the  large  eddy  is  transferred  from  a  detection  of  a  typical 
large  coherent  structure  to  its  production.  The  turbulent  spot  may  be  a 
basic  module  contributing  to  the  oreation  of  a  turbulent  boundary  layer, 
but  its  composition  is  by  no  means  simple,  because  it  contains  a  hierarchy 
of  eddies  arranged  in  a  particular  order,  rather  than  a  single  eddy  which 
scales  with  the  boundary  layer  thickness.  Thus,  the  precise  relationship 
between  the  turbulent  spot  and  the  large  coherent  eddy  remains  to  be  shown, 
although  it  is  known  that  the  spot  retains  some  of  its  features  indefinite¬ 
ly  in  spite  of  the  prolonged  interaction  with  a  turbulent  boundary  layer 
(Zilberman,  Hygnanski  and  Kaplan  1977;  Haritonidis,  Kaplan  and  Wygnanski 
1977;  hereafter  referred  to  as  ZWK  and  HKW,  respectively). 

The  onset  of  turbulence  in  a  boundary  layer  does  not  occur  along  a 
continuous  front  but  rather  at  isolated  spots  which  spread  while  being  con- 
vected  downstream,  eventually  coalescing  with  one  another  to  form  a  turbu- 
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lent  boundary  layer  (Emmons  1951)*  The  spots  occur  randomly  in  time  and 
spaoe  but  they  are  easily  triggered  by  roughness  elements  surface  imperfec¬ 
tions,  small  jets  or  other  disturbances.  In  view  of  their  arrowhead  shape 
the  spots  are  the  most  clearly  identifiable  features  in  the  transition  pro¬ 
cess.  They  may  be  linked  to  the  instability  of  the  laminar  boundary  layer 
(Gaster  1968;  Wygnanski,  Haritonidis  and  Kaplan  1979)  on  one  hand  and  to 
the  turbulent  boundary  layer  on  the  other.  The  purpose  of  this  investiga¬ 
tion  is  to  provide  quantitative  data  describing  the  relationship  between 
interacting  spots  and  fully  turbulent  boundary  layer. 

The  spot  was  discovered  by  Emmons  (1951)  who  also  developed  proba¬ 
bilistic  calculations  for  the  transition  process  based  on  the  following  ob¬ 
servations  and  assumptions: 

(i)  The  spot  originates  at  a  point. 

(ii)  There  is  a  sharp  distinction  between  the  turbulent  re¬ 
gion  within  the  spot  and  in  laminar  surroundings. 

(iii)  The  spot  grows  uniformly  while  being  converted  down¬ 
stream  . 

(iv)  There  is  no  interaction  among  spots. 

On  this  basis  an  expression  was  obtained  by  Emmons,  predicting  the  in- 
termittency  at  any  point  in  the  flow  and  enabling  calculation  of  sxin  fric¬ 
tion  in  the  transition  region. 


Some  features  of  the  turbulent  spot  were  first  documented  by  Schubauer 
and  Klebanoff  (1956).  The  spot  was  initiated  by  means  of  am  electric  dis¬ 
charge  in  a  laminar  boundary  layer  and  was  detected  by  a  hot  wire  located 
in  the  boundary  layer  downstream  of  the  spark.  The  general  shape  of  the 
spot,  the  celerity  of  its  boundaries  and  the  angles  at  which  it  spreads, 
can  all  be  found  in  their  report.  They  also  observed  that  a  fixed  three 
dimensional  roughness  element  located  on  the  surface  produces  a  turbulent 
wedge  in  its  wake  which  spreads  laterally  at  the  same  angle  as  the  spot, 
suggesting  that  the  turbulent  wedge  may  have  been  generated  by  a  succession 
of  spots  telescoping  into  one  another. 

Elder  (I960),  who  used  the  same  technique  to  initiate  the  spots  as 
Scubauer  and  Klebanoff,  provided  an  experimental  evidence  for  the  fourth 
assumption  of  Emmons.  He  simultaneously  generated  two  spots  at  different 
spanwise  locations  and  measured  the  resulting  intermittency  factor  down¬ 
stream,  concluding  that  the  area  covered  by  turbulence  was  simply  the  com¬ 
bined  area  of  the  spots,  assuming  that  they  grew  independently  of  one 
another.  Coles  and  Barker  (1975)  made  their  observations  in  water  using  a 
laser  doppler  velocimeter.  They  concluded  that  the  ensemble  averaged  tur¬ 
bulent  spot  in  a  laminar  boundary  layer  can  be  represented  by  a  single 
large  vortex,  which  entrains  irrotational  fluid  from  the  free  stream  above 
it  and  rotational  fluid  from  laminar  boundary  layer  surrounding  it.  An  at¬ 
tempt  was  made  to  generate  a  synthetic  turbulent  boundary  layer  by  trigger¬ 
ing  several  spots  simultaneously  and  noting  the  resulting  mean  velocity 
profile.  Savas  (1979)  constructed  a  synthetic  turbulent  boundary  layer  in 
air,  by  periodically  disturbing  the  flow  with  protruding  arrays  of  rough- 
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ness  elements,  which  were  placed  at  various  locations  along  the  span. 
Savas  observed  that  some  spots  shifted  their  position  laterally  as  a  result 
of  an  interaction  with  their  neighbours.  The  lateral  "transposition"  of 
spots  might  have  been  associated  with  a  natural  regeneration  of  new  struc¬ 
tures  at  the  wing  tips  of  the  artificially  evoked  spots  (Wygnanski,  Harito- 
nidis  and  Kaplan,  1979;  hereafter  referred  to  as  WHK). 

Wygnanski,  Sokolkov  and  Friecta an  (1976;  hereafter  referred  to  as  WSF) 
mapped  the  flow  field  within  the  spot  by  measuring  all  three  components  of 
the  ensemble  averaged  velocity  inside  it.  They  checked  the  shape  of  the 
spot,  its  spreading  angles  and  the  scales  governing  its  similarity  with 
downstream  distance.  Furthermore,  they  proved  that  the  spot  attains  a  un¬ 
iversal  structure  and  shape  regardless  of  the  type  and  intensity  of  the 
disturbance  generating  it. 

The  prevailing  views  on  the:  "  Coherent  Structure  of  Turbulent  Boun¬ 
dary  Layers"  are  expressed  in  a  volune  of  proceedings  from  a  workshop  dedi¬ 
cated  to  this  subject  in  addition  to  a  number  of  survey  articles  which  can 
be  found  in  the  open  literature  (Laufer  1975,  Willmarth  1975,  Willmarth  and 
Bogar  1977  etc.).  The  views  held  most  commonly  are  backed  by  flow  visuali¬ 
zation  and  by  measurements  using  conditioned  sampling  techniques.  Laufer 
(1975)  claims  that  the  double-layered  nature  of  the  turbulent  boundary 
layer  may  owe  its  existence  to  two  types  of  coherent  structures.  The  outer 
layer  oontains  large  scale  three  dimensional  vorticity  lumps,  while  the 
inner  layer  may  contain  a  thin  vortex  sheet  which  is  occasionally  lifted  by 
the  outer  structures,  resulting  in  the  generation  of  wall  streaks  (see  also 


Blackwelder  1978).  The  relation  between  the  inner  and  outer  structures  was 
Investigated  by  Kin  et  al.(1971),  Blackwelder  (1978)  and  others.  The  lift¬ 
ing  of  the  wall-streaks  occurs  periodically  and  is  often  referred  to  as 
"bursts";  although  the  bursts  occur  near  the  wall  the  bursting  periods 
scale  with  the  outer  flow  parameters  rather  than  with  parameters  governing 
the  wall  region  (Rao  et  al.  1971). 

Visual  studies  (Kline  1978,  and  Gupta  et  al.  1971)  suggested  that 
the  spacing  between  streaks  scales  with  the  wall  coordinates,  in  spite  of 
the  fact  that  the  outer  bulges  are  an  order  of  magnitude  wider  than  the 
streaks.  Faleo  (1977,  78)  subdivided  the  eddies  in  the  outer  and  inner  re¬ 
gions  of  the  turbulent  boundary  layer  into  two  categories.  In  the  outer 
region  he  calls  them  large  eddies  and  "typical  eddies",  while  in  the  inner 
region  they  are  referred  to  as  "pockets"  and  "streaks".  The  different 
names  refer  to  different  shapes  of  eddies  which  presumably  also  contain 
different  levels  of  energy.  The  proliferation  of  names  given  to  visual  ob¬ 
servations  attempting  to  describe  kinematically  a  physical  process  stresses 
our  limited  understanding  of  this  subject. 

The  possibility  that  the  transitional  spot  and  the  large  coherent  eddy 
in  a  turbulent  boundary  layer  may  have  something  in  common  was  explored  by 
ZWK  (1977).  A  spot  was  initiated  among  spherical  roughness  elements  which 
tripped  the  boundary  layer  and  was  allowed  to  interact  with  this  turbulent 
boundary  layer.  The  evoked  structure  was  then  tracked  downstream  up  to  120 
typical  turbulent  boundary  layer  thicknesses  (see  also  HKW  1977).  The 
structure  evoked  in  this  manner  agrees  in  detail  with  the  structure  ob- 
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served  In  the  outer  region  of  the  turbulent  boundary  layer  by  Kovasznay, 
Kibens  and  Blaokwelder  (1970)  and  Justifies  the  effort  invested  in  studying 
turbulent  spots.  The  scale  of  the  artificially  evoked  structure  near  the 
nail  is  approximately  10<$  in  the  streamwise  direction  becoming  2-35  near 
the  interface.  The  spanwise  extent  of  the  evoked  structure  is  less  than 
46  ;  and  its  celerity  is  0.911^.  HKW  (1977)  repeated  the  experiment  using 
a  rake  of  hot-wire  anemometers  and  a  more  powerful  data  acquisition  system. 
They  extended  the  range  of  measurements  and  were  capable  of  assessing  the 
distortion  of  the  evoked  structure  with  downstream  distance.  The  relation¬ 
ship  between  the  educed  eddy  in  a  turbulent  boundary  layer  and  the  artifi¬ 
cially  evoked  spot  in  a  laminar  boundary  layer  is  still  not  well  understood 
in  view  of  the  disparity  in  soales  of  the  two  structures.  Furthermore, 
both  structures  contain  finer  soales  about  whose  origin  and  dynamics  we 
know  rather  little. 


Cantwell,  Coles  and  Dimotakis  (1978;  hereafter  referred  to  as  CCD) 
attribute  the  additional  fragmentation  of  the  spot  (or  outer  bulge)  to  a 
centrifugal  instability  of  the  Taylor-Gortler  type.  Klebanoff,  Tidstrom 
and  Sargent  (1962)  attribute  the  generation  of  small  scale  turbulence  to 
the  breakdown  of  Tollmien  Sohlichting  waves.  WHK  (1979)  observed  a  pair  of 
oblique  wave  packets  located  behind  the  wing-tip  of  the  spot  with  their 
wave  number  vector  being  inclined  at  50°  to  the  streamwise  direction. 
These  wave  packets  were  identified  as  Tollmien-Schliohting  waves  and  their 
breakdown  at  a  higher  Reynolds  number  caused  a  generation  of  new  spots. 
Gad  el  Hak,  Blaokwelder  and  Riley  (1980;  hereafter  referred  to  as  GBR) 
also  suggest  that  the  dominant  mechanism  contributi.ig  to  the  lateral  growth 
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of  a  turbulent  spot  is  caused  by  distabilization  of  the  surrounding  laminar 
boundary  layer. 

In  an  attempt  to  understand  the  process  by  which  transitional  spots 
merge  and  generate  a  turbulent  boundary  layer,  a  controlled  experiment  was 
undertaken  which  simulates  the  interaction  of  spots.  Some  of  the  prelimi¬ 
nary  conclusions  drawn  from  the  interaction  of  successive  spots  originating 
from  a  simple  point  source  were  reported  by  Wygnanski  (1978).  The  evolu¬ 
tion  of  two  spots  evoked  simultaneously  at  different  spanwise  locations  was 
also  investigated,  in  order  to  supplement  the  data  of  Elder  (I960).  A 
preliminary  study  of  spots,  originating  at  different  spanwise  locations  and 
different  times,  was  also  made.  The  broad  aim  of  these  experiments  was  to 
identify  and  compare  the  characteristic  features  observed  in  isolated  spots 
with  the  similar  features  occurring  during  the  interaction  of  spots  and  in 
a  fully  turbulent  boundary  layer.  Hopefully  this  comparison  would  contri¬ 
bute  to  the  understanding  of  this  complicated  flow. 


CHAPTER  B 
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b.i  iaaliuaaai  (fig.  b.i.d. 

All  measurements  and  calibration  procedures  were  conducted  in  the  low 
speed,  low  turbulence,  closed  circuit  wind  tunnel  which  was  constructed  by 
Kenney  Engineering  in  Monrovia,  California.  The  test  section  of  the  tunnel 
is  61  cm  wide  610  cm  long  and  91  cm  high.  The  floor  and  ceiling  of  the 
test  section  are  made  of  plexiglass  and  are  mounted  on  jacks,  thus  enabling 
the  height  of  the  section  to  be  changed  from  90  cm  to  152  cm. 

The  wind  is  generated  by  an  axial  blower  powered  by  a  25  H.P  DC  motor 
and  a  variable  speed  controller.  The  velooity  in  the  test  section  in  ab¬ 
sence  of  any  disturbance  can  vary  from  0  to  a  maximum  of  45  m/see.  The 
maximum  non  uniformity  of  speed  across  the  test  section  is  less  than  0.5< 
exclusive  of  the  boundary  layer  region.  The  longitudinal  turbulence  level 
in  an  empty  test  section  is  approximately  0.04J  at  30  m/sec. 

B.2  The  Plfrte 

An  aluminium  plate,  360  cm  long,  90  cm  wide  and  0.4  cm  thick,  was 
mounted  vertically  in  the  test  section  28  cm  downstream  of  the  inlet  and  11 
cm  from  the  back  tunnel  wall  (fig.  B.2. la).  The  plate  is  pivoted  at  the 
leading  edge  and  can  be  inclined  at  small  angles  of  attack  relative  to  the 
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flow.  Adjustable  screws  at  the  back  of  the  plate  hold  the  plate  in  posi¬ 
tion  relative  to  the  wall  of  the  test  section;  these  screws  can  also  serve 
to  remove  residual  waviness  from  the  surface  of  the  plate. 

The  front  side  upstream  part  of  the  plate  was  machined  to  have  a  cir¬ 
cular  section  of  2  cm  diameter  where  the  front  surface  was  tangent  to  that 
(fig.  B.2.1b).  A  60  cm  adjustable  flap  was  installed  at  the  trailing  edge 
of  the  plate  in  order  to  control  the  circulation  and  assure  that  the  stag¬ 
nation  point  is  located  on  the  test  surface.  A  slight  divergence  of  the 
top  and  bottom  walls  of  the  tunnel  offered  another  degree  of  freedom  for 
fine  adjustment  of  pressure  distribution,  narrow  slots  on  both  sides  of  the 
plate  (Amini  1978)  in  addition  to  the  corner  fillets  were  left  in  order  to 
reduce  the  influence  of  the  corners.  Pressure  measurements  on  the  surface 
of  the  plate  showed  a  maximum  residual  gradient  of  : 

2  -5-1 

(2/pU  )x( dp/dx)  =  -5x10  cm  at  an  ambient  velocity  of  10  m/s.  The  velo¬ 
city  profile  in  the  developing  boundary  layer  was  measured  at  five  differ¬ 
ent  locations  downstream  of  the  leading  edge  and  compared  with  the  theoret¬ 
ical  profile  of  Blasius  (fig.  B.2.2). 

Whenever  a  comparison  with  a  turbulent  boundary  layer  was  required, 
the  boundary  layer  tripped  by  the  same  manner  as  that  of  ZWK  (1977);  by 
rows  of  spheres  1.5  mm  in  diameter  which  were  distributed  in  a  spanwise  di- 
rectlon  at  an  x  location  corresponding  to  Re5  =500.  Turbulent  spots  were 
initiated  by  an  electric  discharge  generated  by  a  commercial  car  ignition 
unit  and  driven  by  a  signal  generator  which  was  controlled  by  the  computer. 
Two  electrodes  0.5  cm  apart  were  inserted  into  a  plexiglass  plug  which  was 
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mounted  flush  with  the  surface;  the  plug  served  as  an  insulator  (fig. 
B.2.1b). 

b.3  Xte  IrAYftr.slnft  Msshaniaa 

The  traversing  mechanism  was  designed  to  translate  the  probes  in  3  di¬ 
rections  x,y,z  and  rotate  them  in  the  X-Y  plane  (fig.  B.3.1).  The  move¬ 
ments  in  the  spanwise  (Z)  direction  and  in  the  streamwise  (X)  direction 
along  the  plate  are  operated  manually  with  an  accuracy  of  0.5  mm.  The 
movements  in  the  Y  direction  (normal  to  the  surface  of  the  plate)  and  the 
rotation  in  the  X-Y  plane  were  controlled  by  the  computer  via  SLO-SYN 
M061-FD08  stepping  motor  and  STM  1800  DV  controller.  The  repeatability  of 
the  motion  in  the  Y  direction  was  approximately  0.04  mm.  The  angular  mo¬ 
tion  covered  a  range  of  45°  with  an  accuracy  and  repeatability  of  1/2°. 
One  cannot  translate  the  probe  while  rotating  it.  The  mode  of  movements  is 
selected  manually  by  tightening  screw  number  1  and  releasing  screw  number  3 
for  translation  or  by  tightening  screw  number  3  and  releasing  screw  number 
1  for  rotation. 

b.4 

Multi  channel  constant  temperature  hot  wire  anemometers  were  con¬ 
structed  in  the  Tel-Aviv  University  electronic  shop  from  a  basic  circuit 
designed  by  R.E.  Kaplan  and  used  at  the  University  of  Southern  California 
(fig.  B.4.1).  Two  sets  of  variable  potentiometers  provided  a  choice  of 


over-heat  ratios  ranging  from  0.05  to  1.  Present  measurements  were  made  at 
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0.5  over-heat  ratio;  the  frequency  response  of  the  set  was  better  than  10 


KHz,  for  5y  diameter  tungsten  wires. 


b.5  flat,  flics,  flaflaa  Iiaaflaal 


Most  of  the  results  presented  in  this  work  were  obtained  using 
home-made  hot  wire  rakes  of  two  kinds  :  (i)  10  normal  wires,  parallel  to 

the  surface  of  the  plate  and  normal  to  the  flow,  (ii)  4  X-wires  inclined 
at  45°  to  the  surface  and  to  the  flow. 


Twenty  jewellers  broaches  were  glued  in  pairs  onto  a  fiberglass  board, 
5  cm  long,  1  cm  wide  and  0.1  cm  thick.  Each  pair  of  broaches  served  as  a 
prong  onto  which  a  5ptungsten  wire  was  welded.  The  distances  of  the  wires 
from  the  bottom  surface  were  1.35,  2.39,  3.42,  4.52,  6.28,  8.63,  11.23, 
15.48,  20.18  and  24.88  cm  (fig.  B.5.1). 


In  the  probe  of  the  second  kind,  sixteen  broaches  were  connected  to 
the  board  in  a  configuration  yielding  a  rake  of  4  X-wires.  Tne  distance 
between  each  slanted  wire  and  its  neighbour  was  about  1  mm  and  between  the 
centers  of  each  X-wire  it  was  about  2  mm  (fig.  B.5.2). 
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CHAPTER  C 

mi  Acauismofl 

The  digitizing  of  the  data  was  done  using  a  minicomputer  (VAR1AN  72 
model)  having  a  32  K  word  memory,  16  bit  each.  Seven  and  nine  track  digi¬ 
tal  tapes  were  available  for  storing  the  data.  The  data  was  digitized  with 
12  bit  precision  (4096  quantizing  steps).  Data  acquired  by  normal  wires 
was  converted  to  velocity  and  recorded  on  digital  tapes  immediately.  Data 
acquired  with  X-wires  was  digitized  and  recorded  in  its  raw  form,  this  data 
was  later  processed  on  a  PDP  11/60  computer  which  had  a  memory  of  128  K 
words,  2  RK07  disks  (27  M  byte  each)  and  was  thus  more  convenient  to  use 
than  the  VARIAN. 

c.  i  ilia  .Qalitoatloa  &£  Normal  licaa 

The  hot  wires  were  calibrated  by  fitting  a  third  order  polynomial  to 
four  calibration  velocities,  in  the  range  1  m/s  to  10  m/s.  The  latter  vel¬ 
ocity  corresponded  to  the  free  stream  velocity  outside  the  boundary  layer. 
Thus  for  a  given  voltage  E^^  (where  the  subscript  i  denotes  the  ith  wire  in 
the  array)  the  velocity  Uj  is  : 

(1)  -  *1«B1x(E1-Eot).C1x(E1-Eol)2.D1x(El-Eol)3 

where  A^,  B^,  C^,  and  D^^  are  the  four  calibration  coefficients  and  is 

the  voltage  measured  in  the  abseenee  of  any  velocity,  i  is  the  index  of 
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the  hot  wire  probe  (i=1 , . . ., 10);  j  is  the  index  of  calibration  velocity 
(jsl . 4). 

At  times  when  the  temperature  in  the  laboratory  remained  constant  the 
hot  wires  did  not  drift  by  more  than  2 %  during  an  entire  day. 

C.2  Look-Up  Table  (high  speed  conversion  of  voltages  to  velocities) 

The  concept  of  a  look-up  table  was  suggested  by  Dr.  J.  Haritonidis 
(U.S.J)  while  discussing  various  possibilities  of  X-wires  calibration. 
Since  the  analog  to  digital  converter  has  only  4096  quantizing  steps,  it  is 
possible  to  translate  each  step  to  its  corresponding  velocity,  once  for 
every  calibration.  Hence,  the  conversion  of  voltages  to  velocities  becomes 
a  simple  readout  of  tables  where  the  voltages  are  actually  the  indices  of 
the  vector  terms.  If  for  an  A/D  converter  the  voltage  range  is  from  -V  to 
+¥  volts,  then  each  quantizing  step  corresponds  to  a  number  between  1  and 
4096  with  the  resolution  of  2V/4096  Volts.  Having  four  coefficients  A,  B, 
C  and  D  for  each  wire,  the  voltages  are  translated  to  4096  velocity  values 
by: 

(2)  U(J)  s  A+Jx( B+Jx( JxD+C) ) 
where  J  goes  from  -2048  to  2048 

Usually,  the  number  of  terms  is  smaller  since  not  all  the  range  is  used. 
Finally,  during  'look-up'  conversion  the  single  operation  needed  is  to 
equate  the  velocity  U  of  voltage  J  with  U(J)  term. 
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c.3  ins  fiallteafclflo  at  fcMlwa 


The  calibration  of  a  rake  of  X-wires  has  to  be  efficient  in  order  to 
save  on  computing  time  and  memory  space.  The  fastest  calibration  procedure 
is  attained  through  the  use  of  look-up  tables  having  two  indices  for  the 
velocities  and  the  yaw  angles  of  the  wires.  However,  if  I  and  J  are  the 
quantizing  values  for  the  voltages  and  of  both  wires,  then  a  two  di- 
mensional  look-up  table  would  require  1.6x10  memory  locations  per  X-wire. 
This  number  is  prohibitive,  bearing  in  mind  that  it  has  to  be  multiplied  by 
the  number  of  X-wires  in  the  rake;  therefore,  if  a  look-up  table  concept 
is  to  be  used,  it  must  be  accompanied  by  a  space  saving  scheme.  The  cali¬ 
bration  procedure  is  based  on  the  fact  that  within  the  calibration  range 
the  velocity  vector  Q  and  its  direction  a  can  be  obtained  by  a  transforma¬ 
tion  of  the  two  independent  X  wire  voltages  E1  and  E 2  (see  Willmarth  and 
Bogar,  1977;  Coles,  Cantwell  and  Wadcock,  1978). 

It  is  well  known  that  the  voltage  E  supplied  to  maintain  the  hot  wire 
at  a  oonstant  temperature  is  increasing  monotonically  with  velocity  and  de¬ 
pends  roughly  on  the  cosine  of  the  angle  between  the  wire  and  the  flow  di¬ 
rection.  The  response  of  the  wire  to  a  change  in  the  flow  direction  is  un¬ 
ique  for  90°  only. 

The  relation  between  the  output  voltage  of  each  wire  and  the  appropri¬ 
ate  velocity  and  direction  can  be  expressed  formally  by  the  following  equa¬ 


tions: 
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(3)  E1  =  Ej(Q,a)  ;  Eg  =  EgfQ.cO 

where  Q  is  the  velocity  vector  and  a  is  its  direction  with  respect  to  the 
probe  axis.  From  the  monotonia  behaviour  of  each  wire  it  is  clear  that 
for  every  combination  of  E^  and  Eg  within  the  calibration  region,  one  ob¬ 
tains  only  one  combination  of  Q  and  a.  In  fact,  the  calibration  region  is 
defined  by  the  inverse  transformation  that  obeys  the  relations: 

(4)  Q  =  QCE^Eg)  ;  a  =a(E^  ,Eg) 

In  order  to  predict  the  calibration  range  in  the  E^Eg  plane,  it  is  helpful 
to  assume  a  monotonic  behaviour  of  the  wire  response  for  three  possible 
cases  of  X-wire  configurations  (illustrated  in  figures  C.3>1a,  b  and  c). 
In  the  first  oase,  the  two  slanted  wires  are  perpendicular  to  one  another. 
Therefore,  a  change  in  direction  at  a  constant  velocity  Q  from  “l  t0  “2 
would  cause  a  decrease  in  Eg  and  a  concomitant  increase  in  E1  (fig. 
C.3.2a).  The  curves  for  a  constant  Q  on  E^-  Eg  plane  will  have  the  shapes 
shown  in  figures  C.3.2b  and  c  (when  the  angle  between  the  wires  is  smaller 
or  larger  than  90°,  respectively).  An  inorease  in  the  velocity  will  shift 
the  curve  diagonally  upwards.  At  the  large  angles  of  yaw  the  response  of 
the  wires  may  not  be  unique  as  a  result  of  strong  end  effects,  resulting 
from  the  finite  length  of  the  wires  and  the  prongs,  when  the  flow  becomes 
parallel  to  one  of  the  wires.  Obviously  one  cannot  rotate  the  wire  rela¬ 
tive  to  the  flow  by  an  angle  larger  than  90°,  thus  reversing  the  flow  di¬ 
rection  relative  to  the  particular  wire. 


The  velooity  Q,  measured  by  an  X-wire,  oan  be  interpolated  in  a  three 
dimensional  spaoe  with  two  independent  variables  E.^  and  Eg  in  the  same 
manner  as  the  calibration  curve  for  a  single  wire  is  interpolated  by  the 
third  order  polinomial  having  a  single  variable  E.  Assuming  that  the  func¬ 
tional  dependence  on  the  angle  is  not  more  complicated  than  on  the  veloci¬ 
ty,  the  calibration  procedure  for  the  X-wire  takes  the  following  steps  : 

1 )  Storing  descrete  points  measured  at  selected  velocities  and  angles  - 
plotting  equation  3. 

2)  Calculating  surfaces  for  Q  and  awhile  finding  equation  4. 

3)  Data  reduction  obtaining  the  velocity  Q  and  the  angle  a  by  using  equation 
number  4. 

Since  the  simple  minded  concept  of  the  'look-up*  table  can  not  be  ap¬ 
plied  directly  to  a  rake  of  X-wires,  the  following  procedure  was  used: 

a)  The  E^E,,  planes  were  subdivided  into  60x60  grid  elements.  This  number 
was  chosen  arbitrarily  although  the  accuracy  is  improved  as  the  number  of 
elements  increases. 

b)  The  velocity  and  angle  at  each  grid  point  (figures  C.3.3a  and  C.3.3b) 
were  calculated  by  using  spline  interpolation. 

o)  The  velooity  range  inside  each  square  was  of  the  order  of  15  em/s  cor¬ 
responding  to  1.531  from  the  free  stream  velocity  used  in  this  experiment. 
Similar  angular  resolution  was  not  considered  satisfactory,  because  the 
perturbation  in  the  normal  component  of  velocity  during  the  passage  of  the 
spot  is  much  smaller  than  in  the  streamwise  component,  requiring  a  better 
angular  resolution.  Therefore,  the  values  within  each  element  were  calcu¬ 
lated  by  interpolation,  using  a  first  order  two  dimensional  polynomial. 
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Those  values  were  stored  in  the  form  of  a  table  to  be  read  out  during  data 
reduction.  If  Q  is  the  velocity  inside  a  grid  element  (i,J)  and  (DE^)^j  is 
the  voltage  increment  for  one  wire,  while  (DE^)^  is  the  voltage  increment 
for  the  other  wire;  the  polynomial  describing  the  velocity  is  given  by  : 


(5)  Q  =  Qij+Aij^DVij^ij^Vij^ij^ViJ^^aJiJ 


similar  results  hold  for  the  angles. 

d)  In  converting  voltages  to  velocities  using  the  ’look-up'  table,  the  vol¬ 
tage  for  each  wire,  divided  by  the  width  of  the  element,  is  in  fact  the 
element  number. 


The  determination  of  velocity  and  angle  required  7  divisions  or  multi¬ 
plications  and  4  subtractions  or  additions.  The  above  procedure  is  sche¬ 
matically  described  in  the  Appendix. 


The  present  X-wire  calibration  method  offers  several  opportunities  for 
checking  the  validity  of  the  calibration  range.  Bearing  in  mind  the  limi¬ 
tations  concerning  the  orientation  of  the  slanted  wires  (fig.  C.3.1),  it 
is  recommended  that  before  each  experiment  some  typical  data  records  would 
be  displayed  on  the  calibration  map.  In  figure  C.3.4  an  example  of  E1  and 
E2  voltages  during  the  passage  of  a  spot  in  a  laminar  boundary  layer  at 
four  locations  is  given,  showing  also  the  borders  of  the  calibration  re¬ 
gion.  Such  a  plot  cannot  guarantee  that  all  the  other  events  will  behave 
precisely  in  the  same  manner  but  it  can  show  the  trend  in  each  kind  of 
measurement.  One  may  also  use  the  computer  to  indicate  whenever  a  given 
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voltage  is  outside  of  the  calibration  window  (see  the  Appendix)  so  that  the 
latter  could  be  extended  whenever  the  calibration  window  is  exceeded  by  an 
arbitrary  prescribed  percentage  of  events;  however,  if  the  calibration 
window  is  exceeded  very  rarely,  the  particular  event  may  simply  be  deleted 
from  the  ensemble-averaged  data. 

c.4  jjgdaaUfln. 

In  all  cases  considered,  data  was  obtained  with  a  rake  of  normal  wires 
as  well  as  with  X  arrays.  The  normal  probe  rake,  which  spanned  the  entire 
boundary  layer,  provided  instantaneous  readings  of  streamwise  velocity  com¬ 
ponent  in  a  single  event  throughout  the  spot,  while  the  multichannel  X 
probe  provided  a  two  dimensional  flow  field  over  a  portion  of  the  spot. 
With  neither  probe  one  was  able  to  study  the  boundary  layer  very  close  to 
the  wall.  The  closest  wire  to  the  wall  for  the  normal  rake  was  0.07  6fc  and 
for  the  X-wire  rake  0.12  6^.  The  spatial  resolution  of  the  X-wire  is  lim¬ 
iting  its  usefulness  and  accuracy  whenever  the  velocity  gradient  is  strong. 
Unfortunately  this  fact  is  inescapable  in  a  normal  study  of  a  boundary 
layer.  At  the  measuring  location  closest  to  the  wall  the  mean  velocity 
gradient  did  not  exceed  1  m/s  per  1  mm  resulting  in  a  possible  error  for 
the  mean  velocity  which  could  not  exceed  3-4$.  The  streamwise  component  of 
velocity  measured  with  an  X  array  was  compared  to  that  measured  with  a  rake 
of  normal  wires,  giving  additional  confidence  to  the  data.  Figure  C.4.1 
shows  the  ensemble-average  of  the  streamwise  perturbation  velocity  compo¬ 
nent  during  a  passage  of  a  transitional  spot  in  a  laminar  boundary  layer. 
The  dotted  lines  represent  the  time  history  of  the  streamwise  component  of 


velocity  obtained  from  the  normal  wires,  and  the  solid  lines  represent  the 
same  data  obtained  from  X-wires.  It  is  clearly  seen  that,  within  the  accu¬ 
racy  mentioned  before,  the  agreement  between  the  two  sets  of  data  is  good 
indeed.  The  distribution  of  u*  measured  in  the  turbulent  boundary  layer  is 
compared  in  figure  C.4.2,  for  data  obtained  by  the  rake  of  X-wires  and  by 
the  rake  of  normal  wires.  At  y/^.=0.12  the  measurements  with  an  X  array 
show  that  the  intensity  of  the  streamwise  velocity  component  is  6.5>  small¬ 
er  than  the  intensity  measured  with  an  array  of  normal  wires.  The  differ¬ 
ence  decreases  further  away  from  the  wall  and  the  intensity  deduced  from 
X-wire  measurements  becomes  higher  in  about  5%  at  y/6fc . 

In  order  to  obtain  statistically  meaningful  results,  each  measurement 
was  repeated  between  100  and  300  times  depending  on  the  type  of  data  de¬ 
sired.  Each  record  contained  770  data  words  per  wire,  which  included  two 
identification  words  at  the  end.  The  sampling  rate  was  4000  per  second  per 
channel  corresponding  to  time  resolution  of  250  ysee  and  the  skew  time 
between  wires  was  about  10  ysee.  Each  record  thus  contained  192  msec  of 
data  at  frequencies  not  exceeding  2  KHz.  In  order  to  represent  the  veloci¬ 
ty  as  a  function  of  time,  a  use  of  a  one  dimensional  array  should  be  suffi¬ 
cient;  however  3  indices  can  describe  the  velocity  at  any  time  at  any  one 
of  the  10  wire  locations  in  any  particular  event.  For  example  the  velocity 
of  wire  number  K  at  time  J  in  event  I  may  be  written  as: 

(6)  U(I, J,K) 

Trying  to  identify  a  general  overall  shape  of  a  repeatable  phenomenon,  one 


Page  21 


may  ensemble-average  the  data  and  express  It  in  the  form 


(7)  <U( J,K)>  =  l  U(I, J,K)/P 

1-1 

P  =  number  of  events. 


Temporally  averaged  data  which  is  only  meaningful  for  steady  flow  may  be 
evaluated  by  sunning  up  the  ensemble-average  for  each  individual  .wire.  By 
doing  so,  one  tacitly  assumes  that  the  duration  is  sufficient  to  give  a  me¬ 
aning  fill  time  averaged  term.  The  number  of  arguments  in  this  case  reduces 
to  one,  i.e  the  wire  number 


(8)  U(K)  =  1  <U( J,K)>/N 

J~1 

N  a  number  of  points. 


However,  when  the  flow  is  unsteady  like  in  the  case  of  an  isolated  spot, 
only  the  ensemble-averaged  of  the  perturbation  velocity  has  a  physical  me¬ 
aning.  Using  the  above  mentioned  formulation,  the  perturbation  velocity  is 
defined  by 


(9)  U  (I,J,K)=U(I,J,K)-U.(K) 
P  1 


where  U^(K)  is  calculated  by  using  (8)  in  the  laminar  boundary  layer.  The 
value  for  the  mean  perturbation  is  obtained  by  combining  (7)and  (9): 


(10)  <UJJ,K)>=  l  Un(I,J,K)/P 


CHAPTER  D 
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Structure 


The  flow  field  in  isolated  transitional  spots  was  remeasured  for  the 
sake  of  comparison  with  the  flow  field  created  by  an  array  of  interacting 
spots  and  by  a  fully  turbulent  boundary  layer.  Furthermore,  new  measuring 
methods  using  a  rake  of  X-  wires  offered  the  opportunity  to  map  the  two 
components  of  the  velocity  vector  across  the  entire  boundary  layer  simul¬ 
taneously.  One  would  like  to  compare  the  single  spot  growing  in  a  laminar 
boundary  layer  and  the  large  coherent  structure  educed  in  a  fully  turbulent 
boundary  layer  after  a  prolonged  interaction  between  the  turbulence  in  a 
transitional  spot  and  the  turbulence  in  the  surrounding  boundary  layer. 
ZWK  (1977)  observed  that  only  a  small  portion  of  the  transitional  spot  re¬ 
tained  its  coherence  as  a  result  of  the  above-mentioned  interaction.  The 
scales  of  the  educed  structure  were  approximately  106  in  the  streamwise  di¬ 
rection  and  about 46  in  the  spanwise  direction.  The  length  and  width  of  the 
transitional  spot  at  a  comparable  location  are  approximately  five  times 
larger. 


Figure  D. 1 . 1  describes  the  mean  perturbation  contours  of  the  stream- 
wise  velocity  component  measured  by  a  rake  of  hot-wires  which  are  normal  to 
the  direction  of  the  mean  flow,  and  located  at  various  y  distances  from  the 
surface.  The  ordinate  in  this  figure  is  stretched  in  order  to  exhibit  the 
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flow  field.  A  proper  height  to  length  ratio  is  obtained  by  multiplying  the 
absoissa  by  a  factor  of  65  after  converting  the  time  axis  to  distance  by 
using  the  free  stream  velocity;  this  number  would  be  reduced  to  52  if  one 
were  to  use  the  convection  velooity  of  the  spot  for  the  same  purpose. 

The  flow  field  described  by  velooity  perturbation  relative  to  the  lam¬ 
inar  boundary  layer  in  the  plane  of  symmetry  of  the  spot  can  be  divided 
into  two  main  regions.  The  rounded  triangular  region  removed  from  the  wall 
has  a  velooity  deficit  in  comparison  with  the  unperturbed  laminar  profile, 
while  the  velocity  close  to  the  surfaoe  is  increased  relative  to  the  unper¬ 
turbed  Blaslus  velocity  profile.  In  both  regions,  the  maximum  values  of 
the  perturbation  velooity  is  about  30>(J  .  A  similar  description  was  alrea- 

oo 

dy  reported  by  ZWK  (1977).  One  may  infer  from  the  literature  (Coles  and 
Barker  197*1,  WSF  1976)  that  the  ensemble-averaged  data  gave  rise  to  the  in¬ 
terpretation  that  the  spot  can  be  regarded  as  a  single  eddy.  Hot  all  the 
features  of  the  spot  can  be  attributed  to  such  a  structure.  Furthermore, 
recent  data  Indicates  that  the  spot  contains  an  orderly  array  of  smaller 
eddies  (Leonard  1979,  Wygnanski  1980);  nevertheless  the  ensemble-averaged 
data  is  helpful  for  modeling  the  flow.  Figure  D.1.2a  Is  a  superposition  of 
the  mean  two-dimensional  velocity  perturbation  vectors  on  the  streamwise 
-1>  and  - 2 i  perturbation  contours  [i.e.  (U-U1)/Uoos-0.01  and  -0.02] 
shown  in  figure  0.1.1.  The  mean  values  of  the  velooity  components  prevail¬ 
ing  in  laminar  flow  in  the  x  and  y  directions  respectively  were  subtracted 
from  the  averaged  velocity  measurements  for  each  measuring  station. 
Hamely,  each  streamwise  component  shown  is  given  by  U-U^ ,  while  each  normal 
component  is  given  by  (V-V.)xC;  U.  and  V,  are  the  mean  laminar  values  of 
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the  two  components.  In  this  way  the  flow  field  before  the  arrival  of  the 
spot  (T  <  10  milliseconds)  is  described  by  dots  which  represent  an  unper¬ 
turbed  region.  The  scale  factor  C  corresponds  to  the  coordinate  stretching 
used  in  figure  D- 1.1  and  facilitates  the  observation  of  the  circulatory  mo¬ 
tion  of  the  flow  field  in  the  spot.  Along  most  of  the  leading  edge  of  the 
spot  (20  <  T  <  30  milliseconds,  and  y  >  4  mm)  the  flow  is  directed  towards 
the  surface,  while  the  perturbation  velocity  of  the  streamwise  component 
becomes  negative.  The  perturbation  velocity  of  the  normal  component  decre¬ 
ases  and  changes  its  sign  near  the  surface  (y=2.5  mm)  at  the  leading  edge 
of  the  spot.  The  perturbation  velocity  of  the  streamwise  component  attains 
a  minimum  in  the  central  region  of  the  spot  (corresponding  to  T=35  msec  and 
y  >  2  nm)  while  the  normal  component  changes  direction.  This  region  may  be 
considered  as  the  core  of  the  large  eddy.  At  a  later  time  (T  >  35  msec) 
the  gradient  of  the  streamwise  velocity  component  dU/dT  becomes  positive. 
The  velocity-defect  region  terminates  at  the  trailing  interface,  but  the 
positive  perturbation  in  V  lingers  on.  The  negative  V  perturbation  compo¬ 
nent  transports  momentum  to  the  fluid  near  the  surface  causing  a  positive 
perturbation  in  the  streamwise  component  of  velocity,  while  a  positive  V 
perturbation  occurring  at  a  later  time  is  associated  with  a  deficit  in  the 
U  component  of  velocity  at  the  outer  region. 

At  the  first  measuring  level  (y=  2.5  mm)  near  the  leading  interface 
(TsIO  msec),  one  may  observe  a  positive  perturbation  in  both  U  and  V  simul¬ 
taneously.  The  positive  perturbation  of  the  normal  velocity  component  (V) 
at  the  lowest  measuring  level  (y=2.5  ram)  is  attributed  to  a  vortex-like  mo¬ 
tion  near  the  front  portion  of  the  spot  ( 15  <  T  <  70  msec),  caused  by  fluid 


arriving  at  the  plane  of  symmetry  from  other  spanwise  locations. 

The  normal  velocity  component  attains  a  maximum  near  the  trailing  edge 
of  the  spot  and  decreases  smoothly  to  its  non-dlsturbed  level  behind  the 
spot.  Very  close  to  the  wall,  the  perturbation  in  the  fluid  motion  is  di¬ 
rected  towards  the  spot's  center.  This  region  which  was  named  by  Schubauer 
and  Klebanoff  (1956)  as  the  recovery  trail,  has  been  shown  to  be  stable  and 
devoid  of  velocity  fluctuations.  Assuming  that  the  perturbation  flow  field 
in  the  ensemble  averaged  spot  results  from  the  presence  of  a  horseshoe  vor¬ 
tex,  the  direction  of  flow  in  the  plane  of  symmetry  near  the  trailing  in¬ 
terface  could  be  determined  by  the  outward  motion  caused  by  a  pair  of 
counter-rotating  strearawise  vortices  which  are  at  the  tail  of  the  horseshoe 
vortex  .  The  outward  motion  results  in  an  additional  region  of  velocity  de¬ 
fect  which  is  represented  by  the  -1*  0  contour  from  figure  D. 1.1.  This 
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region  was  explored  lately  by  Van  Atta ,  Sokolov,  Antonia  and  Chambers 
(1981)  ,  while  they  were  measuring  the  flow  disturbance  produced  in  the 

free  stream  above  the  turbulent  spot . 

The  advantages  of  measurement  with  a  rake  of  X-wires  become  apparent 
by  comparing  figure  D. 1.2a  with  figure  D. 1.2b.  In  the  second  figure  the 
normal  velocity  component  was  calculated  from  the  continuity  equation  as¬ 
suming  that  the  flow  is  two-dimensional  (3W/3zsO),  The  scale  on  both  fig¬ 
ures  is  identical .  The  calculated  V  shows  that  the  overall  direction  of 
flow  within  the  spot  is  similar  to  the  observed  one,  but  not  identical. 
The  angular  motion  in  the  ensemble  averaged  spot  in  both  figures  may  be  mo¬ 
deled  by  a  large  vortex  but  the  detailed  comparison  casts  some  doubt  about 
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the  assumption  of  two  dimensionality.  A  case  in  point  is  the  simultaneous 
positive  velocity  perturbation  of  the  normal  and  streamwise  components 
which  appears  near  the  solid  surface  at  both  interfaces  of  the  spot,  and 
probably  results  from  a  non  vanishing  gradient  term  3W/3z.  Thus,  the  dis¬ 
cussion  of  the  flow  field  during  the  passage  of  the  spot  will  rely  on  the 
measured  V  component  of  velocity,  rather  than  the  two-dimensional, 
ensemble-averaged  continuity  equation. 

d.2  Sn&ainagnS. 

The  spot  entrains  non-turbulent  fluid  from  the  laminar  boundary  layer 
and  from  the  potential  flow  above  it.  There  are  a  number  of  factors  con¬ 
tributing  to  the  entrainment  of  non  turbulent  fluid:  firstly,  there  are  a 
large  number  of  small  eddies  at  the  interface  which  increase  the  contact 
area  between  the  turbulent  and  non  turbulent  fluid  thus  increasing  the  rate 
of  diffusion  of  turbulence;  secondly,  the  non  turbulent  fluid  is  engulfed 
by  the  large  eddies  in  a  similar  manner  as  in  the  free  shear  layer; 
thirdly,  the  large  turbulent  eddies  induce  velocity  perturbations  outside 
the  interface  which  may  further  destabilize  the  laminar  boundary  layer. 
Figure  D.2.1  shows  a  picture  of  a  dyed  spot  photographed  by  GBR  (1980)  in  a 
water  facility.  In  this  picture  one  sees  an  elevation  view  of  eddies  on 
the  plane  of  symmetry  of  the  spot.  One  may  also  distinguish  the  large  par¬ 
cels  of  irrotational  fluid  which  were  engulfed  into  the  spot  between  two 
adjacent  turbulent  eddies.  Similar  observations  were  made  by  Falco  (1977) 
on  a  fully  developed  turbulent  boundary  layer  visualized  by  means  of  smoke 


in  air. 
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In  attempting  to  understand  the  mechanism  which  is  responsible  for  the 
growth  of  the  spot  in  both  spanwise  and  streamwise  directions,  one  has 
first  to  determine  the  region  which  most  actively  entrains  non-turbulent 
fluid.  The  most  obvious  first  step  is  to  consider  the  entrainment  into  an 
ensemble-averaged  spot  on  the  plane  of  symmetry.  For  this  purpose,  one 
should  know  the  shape  of  the  interface  and  the  velocity  field  relative  to 
it.  The  task  is  simple  in  principle  but  requires  tedious  and  difficult 
measurements.  The  first  attempt  was  made  by  WSF  (1976)  who  measured  the 
streamwise  component  of  the  celerity  of  the  leading  and  trailing  interfaces 
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locally  and  calculated  the  function  (U-Uinfc)dy  r#lative  the  respec¬ 
tive  interface.  Inferring  the  entrainment  from  this  procedure  implies  that 
the  flow  relative  to  the  interface  is  steady  and  two  dimensional  [i.e. 
(3W/3z)  a  0].  Furthermore,  the  growth  of  the  spot  in  the  y  direction  has 
been  neglected.  CCD  (1978)  argued  that  in  the  laboratory  coordinates  the 
stream  function  depends  on  x,  y  and  t  and  thus  has  a  total  differential: 
di|»  a  -Vdx  +Udy  +  (3i|>/3t)dt  ;  which  includes  a  time  dependent  term. 
Consequently  the  integration  of  (U-Uinfc)  with  respect  to  y  is  insufficient 
to  calculate  entrainment.  They  estimated  the  rate  of  entrainment  by  assum¬ 
ing  conical  similarity  and  calculating  the  particle  trajectories  relative 
to  the  interface. 


In  the  present  context  an  attempt  is  made  to  calculate  the  rate  of  en¬ 
trainment  directly  from  measurements,  without  resorting  to  the  assumption 
of  two  dimensionality,  the  invariance  of  leading  edge  celerity  with  y,  or 
conical  similarity.  The  average  velocity  field  within  the  spot  and  its  vi¬ 
cinity  was  measured  with  a  rake  of  X-wires.  The  oelerity  of  the  interface 


was  determined  from  the  perturbation  velocity  contours,  by  choosing  ini¬ 
tially  four  characteristic  features  on  the  boundary  of  the  spot  and  deter¬ 
mining  the  time  of  their  occurrence.  These  points  mark  the  following  loca¬ 
tions  : 

locus  1.  The  leading  edge  of  the  spot  near  the  surfaoe. 

locus  2.  The  most  forward-reaching  position  of  the  leading  interface  (the 
overhang).  Both  loci  1  and  2  correspond  to  locus  A  of  CCD. 
locus  3.  The  location  of  the  maximum  height  of  the  spot  (locus  C  of  CCD), 
locus  4.  The  trailing  edge  of  the  spot  near  the  surface  (locus  F  of  CCD). 

By  repeating  the  measurements  at  numerous  streamwise  stations  on  the 
plane  of  symmetry  of  the  spot,  the  loci  labelled  1  through  4  occur  at  dif¬ 
ferent  times  depending  on  the  location  of  the  probe.  These  loci  are  plot¬ 
ted  on  x,t  coordinates,  (fig.  D.2.2),  in  a  manner  suggested  by  CCD.  Data 
from  8  measuring  stations  is  shown  in  figure  D.2.2  for  streamwise  locations 
varying  from  800  to  1500  mm  downstream  of  the  spark.  The  loci  marked  may 
be  connected  by  straight  lines,  whose  slopes  represent  the  streamwise  com¬ 
ponent  of  the  celerity  of  the  particular  features  ohosen.  The  celerity  of 

loci  1  and  2  is  0.9U  and  that  of  locus  3  is  not  appreciably  different 

00 

(0.82U  for  the  Reynolds  number  under  consideration).  The  celerity  of  locus 
00 

4  is  0.55U  .  Loci  1,  2  and  3  are  situated  on  the  leading  interface  of  the 

CO 

spot.  Thus,  one  may  assume  that  the  leading  interface  on  the  plane  of  sym¬ 
metry  moves  downstream  with  approximately  constant  velocity  which  is  inde¬ 
pendent  of  y  (y  being  the  distance  from  the  surface).  Furthermore,  the  in¬ 
clination  and  shape  of  the  leading  interface  does  not  vary  with  increasing 
distance  from  the  spark. 
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The  contours  U-U^/U^s-0.02  which  may  represent  the  boundaries  of  the 
spot  outside  the  laminar  boundary  layer  were  measured  at  U.S.C  by  Wygnanski 
and  Haritonidis  and  are  shown  in  figure  D.2.3*  The  streamwise  separation 
between  adjacent  contours  is  100  mm.  This  determination  of  the  spot’s 
boundary  was  first  used  by  Coles  and  Barker  (1975)  and  is  repeated  here  for 
the  sake  of  convenience;  a  detailed  comparison  between  this  criterion  and 
the  actual  determination  of  the  turbulent-non-turbulent  interface  is  dis¬ 
cussed  by  WSF.  Even  if  the  contour  U-U^/U^,  s-.02  does  not  coincide  with 
the  interface  of  the  spot  it  is  an  impartially  determined  characteristic 
feature  of  the  spot  near  its  outer  boundary.  One  may  deduce  from  figure 
D.2.3  that  the  boundaries  of  adjacent  spots  are  approximately  parallel 
which  could  lead  to  the  notion  that  the  celerity  of  the  trailing  interface 
in  the  X  direction  is  locally  independent  of  y.  In  order  to  compare  the 
celerity  of  the  trailing  interface  over  long  downstream  distances,  one  has 
to  account  for  the  growth  of  the  spot  in  the  Y  direction  as  well.  It 
should  be  remembered,  however,  that  the  coordinates  used  in  figure  D.2.3 
accentuate  the  height  of  the  spot  which  is  physically  a  very  flat  struc¬ 
ture  . 


The  variation  of  the  maximum  height  of  the  spot  with  x  can  be  roughly 
estimated  from  figure  D.2.3;  from  it  one  may  calculate  the  outward  propa¬ 
gation  velocity  of  the  tip  (fig.  0.2.#).  The  forward  'overhang'  of  the 
spot  (locus  2)  does  not  propagate  outwards  in  the  range  of  the  measurements 
shown;  and  it  is  roughly  located  at  y=7  mm  from  the  surface  of  the  plate. 
It  may  be  assumed  that  the  normal  component  of  the  interface-celerity  in¬ 
creases  linearly  with  the  normal  distance  from  the  overhang. 
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Thus: 


(11)  (»)lnt/(*)tlp=K»lnt-«)0,.rluing3/[H-<Y>0,erhang] 


where  H  is  the  maximum  height  of  the  spot. 


Taking  advantage  of  the  general  similarity  of  the  spot  boundaries,  one 
may  calculate  the  streamwise  component  of  the  celerity  of  the  trailing  in¬ 
terface  at  similar  locations  on  the  boundary,  thus  accounting  for  the 
growth  of  the  spot  in  Y  direction.  The  data  shown  in  figure  0.2.3  leads  to 
the  conclusion  that  U.  (trailing  edge)  accelerates  with  increasing  y.  In 
fact,  only  the  lower  portion  of  the  trailing  interface  is  conveeted  at  a 
constant  celerity.  These  results  differ  ffom  the  conclusion  of  WSF  who 
suggested  that  the  oelerity  of  each  interface  is  constant  in  the  plane  of 
symmetry  of  the  spot.  The  difference  may  stem  from  the  fact  that  the  meas¬ 
urements  were  made  locally,  mostly  at  the  lower  part  (y/H  <  0.5)  of  the 
spot  (see  fig.  10  of  WSF). 
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The  celerity  of  the  interface  for  the  elevations  at  which  X-wire  measure¬ 
ments  are  available  are  tabulated  below: 


LEVEL 

yim] 

U.  /U 
le  00 

Table  1 

“te^ 

(¥int/U«>x'° 

1 

2.46 

.9 

.55 

0 

2 

4.31* 

.9 

.55 

0 

3 

6.21 

.9 

.55 

0 

4 

8.21 

.9 

.59 

1.18 

5 

9.96 

.85 

.63 

2.90 

6 

11.84 

.84 

.66 

4.73 

7 

13.71 

.83 

.69 

6.56 

8 

15.71 

.82 

.72 

8.51 

Figure  D.2.5a  represents  the  velocity  vectors  relative  to  the  leading 
interface.  Each  tick  mark  on  the  abscissa  of  this  figure  represents  either 
a  time  scale  equivalent  to  2  msec  or  a  velooity  scale  of  2  m/sec.  Each 
tickmark  on  the  ordinate  represents  a  vertical  distance  equivalent  to  1  mm 
or  a  component  of  velocity  in  m/see  equal  to  1/(U^^x2).  Thus  the  inclina¬ 
tion  of  the  velocity  vector  is  stretched  by  the  same  ratio  as  the  slope  of 
the  interface.  The  vertical  distance  in  figure  D.2.5.a  is  roughly 
stretched  by  a  factor  of  18:1  relative  to  the  horizontal  distance.  The 
same  procedure  was  used  in  figure  D.2.5b  relative  to  the  trailing  inter¬ 
face.  One  may  integrate  the  relative  velocity  component  normal  to  the  in¬ 
terface  and  obtain  the  average  entrainment  distribution  into  the  spot  on 
its  plane  of  symmetry. 
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Table  2 

Measured  Calculated 


Boundary 

Segment 

Segment 

Segment 

Segment 

Segment 

Segmei 

Definition 

1-2 

2-3 

3-4 

1-2 

2-3 

3-4 

Interface 

21} 

19} 

60} 

15} 

5} 

80} 

2}  perturbation 

22} 

14} 

64} 

16} 

6} 

78} 

3}  perturbation 

22} 

22} 

56} 

— 

— 

— 

CCD 

— — i 

... 

... 

16} 

5} 

79} 

The  first  three  rows  in  Table  2  refer  to  the  present  results,  while 
the  results  of  CCD  s>re  shown  for  comparison  in  the  fourth  row.  The  meas¬ 
ured  data  indicates  that  the  leading  interface  entrains  approximately  40} 
of  the  total  entrainment  estimated  this  way.  The  amount  of  fluid  entrained 
under  the  ’overhang’  (segment  1-2)  is  approximately  equal  to  the  amount  of 
fluid  entrained  in  segment  2-3.  The  ratio  of  entrainment  between  the  lead¬ 
ing  and  trailing  interface  corresponds  approximately  to  the  ratio  of  the 
lengths  of  these  ensemble  averaged  interfaces,  suggesting  perhaps  that  both 
are  equally  active  in  the  entrainment  process. 

The  ensemble-averaged  perturbation  velocity  component  normal  to  the 
surface,  which  was  calculated  from  the  U  component  (see  figure  D.1.2b)  by 
assuming  that  the  average  flow  is  two  dimensional  (i.e.,  9W/9z=0),  was  used 
in  order  to  calculate  the  velooity  veotors  relative  to  the  leading  and  tra¬ 
iling  interfaces  which  are  shown  in  figures  D.2.5c  and  d  respectively.  The 
method  of  stretching  the  scale  in  these  figures  is  identical  to  that  used 
in  the  measured  data  (figs.  D.2.5a  and  b)  providing  a  direct  comparison 
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between  the  measured  and  calculated  data.  The  relative  entrainment  by  the 
trailing  Interface  is  much  stronger  if  one  considers  the  calculated  V  rath¬ 
er  than  the  measured  one.  The  calculated  results  are  also  tabulated  in 
Table  2.  The  calculated  distribution  of  entrainment  agrees  fairly  well 
with  the  calculated  results  of  CCD.  The  apparent  inactivity  of  the  upper 
portion  of  the  leading  interface  stems  probably  from  the  fact  that  the  ef¬ 
fects  of  the  surrounding  laminar  boundary  layer  are  neglected  in  the  calcu¬ 
lations. 

D.3  teaaUa  fiaasA  an  a  Steals.  flgaltea&tea 

The  results  which  have  been  shown  so  far  were  analyzed  after  ensemble 
averaging  the  data.  A  hot  wire  rake  containing  nine  wires  made  it  possible 
to  examine  the  time  history  of  the  instantaneous  streamwise  velocity  compo¬ 
nent  across  the  entire  plane  of  symmetry  of  the  spot.  The  velocities  shown 
in  figure  D.3.1  are  expressed  by  equation  9.  The  data  is  filtered  at  250 
Hz  before  the  perturbation  contours  were  calculated  (fig.  D.3. 2).  A  com¬ 
parison  of  the  contours  shown  in  the  present  figure  with  the  corresponding 
ensemble-averaged  contours  shown  in  figure  D. 1.1  is  now  possible.  The  bro¬ 
ken  line  in  figure  D.3. 2  represents  the  +  3$  perturbation  contours. 
Although  the  overall  shape  and  the  instantaneous  contours  are  similar  to 
the  ensemble-averaged  data,  the  internal  structure  appearing  in  a  single 
realization  is  blotted  by  the  averaging  process  for  the  realization  shown 
in  figure  D.3. 2.  One  may  detect,  at  least,  four  regions  (marked  by  'E's  in 
the  figure)  at  which  the  perturbation  levels  attain  a  local  maximum  (mini¬ 
mum).  The  most  negative  perturbation (Up=-0. 35Uto)  appears  near  the  leading 
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edge  of  the  spot  at  an  elevation  corresponding  to  the  local  height  of  the 
laminar  boundary  layer  .  This  region  is  followed  in  time  by  three  addi¬ 
tional  defect  regions  with  perturbation  velocities  ranging  from  -0.3U  to 
-0.1U  . 

CO 

Closer  to  the  surface  there  are  regions  at  which  the  excess  of  veloci¬ 
ty  perturbation  attains  a  local  maximum,  they  usually  occur  at  later  times 
than  the  corresponding  local  minima.  Each  of  these  pairs  of  defect  and  ex¬ 
cess  velocity  regions  may  represent  a  large  coherent  structure.  Present 
observations  indicate  that  the  spot  contains  several  large  eddies  extending 
from  its  outer  boundaries  to  the  surface.  Individual  eddies  were  observed 
previously  in  an  insipient  spot  by  Amini  (1978)  who  also  noted  that  the 
number  of  eddies  increases  with  downstream  distance.  Four  photographs  of 
the  streamwise  velocity -history  taken  by  Amini  (1978)  at  distances  500  >,  x 
>,  200  mm  from  the  perturbation  are  shown  on  figure  D.3.3.  The  free  stream 
velocity  in  this  experiment  was  6  m/s.  The  number  of  eddies  which  may  be 
deduced  from  these  velocity  histories  appears  to  increase  with  x  from  three 
at  x=200  mm  to  five  at  x=500  mm. 


D.4  A  Comparison  Between  Ida  Transitional  Saftt  Aol  Ida  Sacbalaat. 

Baanflacg  Laver  iflgaa  aalaaUil 

The  velocities  measured  in  conjunction  with  isolated  transitional 
spots  evolving  in  a  laminar  boundary  layer  will  now  be  compared  with  con¬ 
ventional  measurements  (time  averaged)  of  velocity  in  a  turbulent  boundary 
layer.  A  sequence  of  ensemble-averaged  velocity  profiles  taken  on  the 
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plane  of  symmetry  of  the  spot  Is  shown  in  figure  0.4.1;  the  mean  profiles 
which  were  measured  at  8  millisecond  intervals  are  represented  by  solid 
lines,  while  the  laminar  profiles  measured  at  the  same  x  location  are 
marked  by  dotted  lines.  The  velocity  perturbation  contours  are  shown  in 
the  lower  part  of  the  figure  for  the  sake  of  comparison.  The  strongest 
velocity  deficit  in  the  profile  oocurs  in  the  central  part  of  the  spot, 
while  the  maximum  velocity  gradient  near  the  wall  and  the  maximum  positive 
velocity  perturbation  occur  at  the  trailing  edge  of  the  spot .  One  may  plot 
the  velocity  profiles  in  the  spot  on  a  logarithmic  scale  and  compare  them 
to  the  universal  profiles  taken  in  a  turbulent  boundary  layer.  Such  a  com¬ 
parison  was  made  previously  by  Coles  and  Barker  (1975),  WSF  (1976)  and  CCD 
(1978).  Six  velocity  profiles  measured  at  X  =103  cm  and  labelled  in  figure 
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D.4.1  are  shown  in  figure  D.4.2.  These  profiles  span  the  entire  duration 
of  the  spot •  A  straight  line  drawn  through  the  data  points  in  the  wall  re¬ 
gion  and  extrapolated  to  y+=10  yields  the  friotion  velocity  U^for  this 
case.  The  same  procedure  was  repeated  at  two  additional  measuring  stations 
84.5  cm  and  122  cm  downstream  of  the  spark  and  the  values  of  U  /U  are  ta- 
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bulated  in  Table  3.  The  time  intervals  between  adjacent  profiles  was  in¬ 
creased  with  increasing  downstream  distanoe  in  order  to  preserve  the  scal¬ 
ing  relative  to  the  overall  duration  of  the  spot .  The  dimensionless  time 
intervals  are  tabulated  in  the  first  row  of  table  3  followed  by  the  values 

.0  /U  for  the  six  profiles  appearing  in  figure  D.4.2.  The  right  hand  column 
x  «*> 

in  table  3  gives  the  friotion  velocities  measured  in  the  tripped  turbulent 
boundary  layer  at  the  same  three  measuring  stations.  The  wake  component 
(Coles  1964)  for  the  profiles  measured  at  X  =103  cm  appear  in  the  fourth 

S 

row  of  Table  3,  and  is  compared  with  the  appropriate  number  in  a  fully  tur- 


bulent  boundary  layer 


Table  3 

T/SPOT  DURATION  0.275  0.375  0.475  0.575  0.675  0.775  T.B.L 

#Xs 


84 

.0385 

.0407 

.0425 

.0445 

.0460 

.0465  .0435 

103 

.037 

.0385 

.0400 

.0415 

.0435 

.045  .  0415 

122 

.037 

.0385 

.0395 

.0410 

.0425 

.0445  .0410 

AU/U  £Xsrl03 

4.8 

4.9 

4.3 

3.9 

3.1 

2.7  2.7 

The  maxim  lira  deviation  of  the  looal  velocity  in  the  outer  part  of  the  turbu¬ 
lent  boundary  layer  is  referred  to  by  Coles  (1964)  as  the  strength  of  the 
wake  component.  This  quantity  (AU/U^)  is  a  function  of  Ree(l.e.  Reynolds 
number  based  on  momentum  thickness).  For  Re0>  6000  the  wake  component 
AU/Ut  equals  2.7  and  it  deoreases  at  lower  Re  until  it  completely  disap¬ 
pears  in  the  neighbourhood  of  Re0«5OO.  When  the  turbulent  boundary  layer 
was  artificially  thickened,  the  strength  of  the  wake,  closer  to  the  distur¬ 
bance,  was  higher  than  its  asymptotic  value  far  downstream.  Coles  attri¬ 
buted  this  effect  to  the  large  eddies  which  originated  at  the  roughness 
elements.  The  data  on  the  last  row  of  Table  3  indicates  that  the  strength 
of  the  wake  component  decreases  towards  the  trailing  edge  of  the  spot,  sug¬ 
gesting  perhaps  that  the  scale  of  the  large  eddies  also  reduces  as  one 
proceeds  towards  the  trailing  edge.  This  observation  is  eonsistant  with 
the  dye  ploture  D.2.1  and  with  the  velooity  contours  shown  in  figure  D.3.2, 
in  which  the  contours  resulting  from  a  single  realization  are  plotted. 
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The  friction  veloeity,  on  the  other  hand,  increases  towards  the  trail¬ 
ing  edge,  reaching  a  maximum  value  at  the  location  corresponding  to  the  en¬ 
semble-averaged  trailing  interface  near  the  surface.  The  friction  veloeity 
in  the  turbulent  boundary  layer  measured  at  the  same  x  location  (X  =103  cm) 

9 

equals  the  friction  velocity  observed  around  the  center  of  the  projection 
of  the  spot  on  the  surface,  and  lags  behind  the  location  at  which  the  spot 
attains  its  maximum  height.  The  average  height  of  the  spot  at  the  location 
at  which  U^is  maximum  is  0.8  H  corresponding  to  1.0  6fc.  The  magnitude  of 
the  wake  component  which  is  equal  to  the  magnitude  prevailing  in  a  turbu¬ 
lent  boundary  layer  occurs  even  closer  to  the  trailing  interface  of  the 
spot.  Identical  observations  were  made  at  other  measuring  stations. 

d.5  felasiu  ElufltoiaLLoaa 

The  measurement  of  veloeity  fluctuations  in  a  non-stationary  flow 
field  requires  careful  definition.  In  a  stationary  turbulent  flow  the  vel¬ 
oeity  fluctuations  are  defined  by: 

(12)  u(J)=  U(J)-U 

where  U(J)  is  the  instantaneous  velocity  at  a  point  J  and: 

T 

(13)  U=  (1/T)|  U(t)dt  . 

0 

In  a  stationary  turbulent  flow  .»e  average  velocity  may  also  be  given  by: 
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(14)  U  =  l  U(J)/N 

J=1 

where  N  is  the  number  of  points  sampled. 

In  the  present  case  the  sampling  of  the  data  was  initiated  prior  to  the  ar¬ 
rival  of  the  spot  and  the  duration  of  sampling  depended  on  the  length  of 
the  spot.  Thus,  the  average  velocity  U  which  depends  on  the  length  of  the 
sampling  window  is  no  longer  an  objective  quantity  and  can  not  be  used  for 
reference.  It  is  customary  to  define  u  by  decomposing  U  into  three  compo¬ 
nents  : 

(15)  U  =  U,  +  <0  >  +  u 

1  P 

where  u  is  the  fluctuating  term;  U  is  the  same  as  in  (9)  and  <U  >  is  the 

1  N  p 

value  of  the  ensemble-averaged  velocity  (<Up>=1/N£  (U-U^) ;N=number  of  real¬ 
izations)  . 

One  may  define  a  local  intensity  of  the  velocity  fluctuations  on  its  root 
ensemble  mean  square  value  (r.e.m.s)  of  the  u  term  (from  eq.  15)  at  each 
measuring  location  as  done  by  Van  Atta  and  Helland  (1980): 

N 

(16)  <u'(J)>  =  [  I  <U(I,J)-U.(I,J)-<U  (I,J)>)2/N]1/2 

1-1  1  p 

where  N  is  the  number  of  events  and  is  calculated  in  the  non-disturbed 


laminar  boundary  layer. 


Page  39 


The  results  computed  from  equation  16  are  shown  in  figure  D.5.1.  Equation 
16  relates  turbulent  activity  to  specific  regions  in  the  spot,  however,  the 
interpretation  of  the  results  described  by  this  equation  is  difficult  and 
prone  to  errors.  The  results  near  the  solid  surface  indicate  that  the  tur¬ 
bulent  intensity  in  the  vicinity  of  the  leading  or  the  trailing  interface 
is  higher  than  in  the  central  region  of  the  sp  t  (Van  Atta  and  Helland 
1980). 

Since  the  length  of  the  spot  is  not  identical  in  every  realization, 
nor  is  the  time  of  arrival  of  the  spot  at  the  measuring  station  precisely 
repeatable,  the  jitter  in  both  quantities  may  contribute  to  an  apparent 
turbulent  intensity.  In  order  to  estimate  this  "apparent"  intensity,  an 
attempt  was  made  to  reduce  the  Jitter  by  aligning  all  events  at  one  of  the 
interfaces.  The  bottom  trace  in  figure  D.5.2  represents  simple 
ensemble-averaged  velocity  in  the  vicinity  of  the  leading  interface.  The 
trace  above  it  was  generated  from  the  same  data;  however,  the  velocity 
signal  in  each  realization  was  shifted  by  the  time  difference  between  the 
arrival  of  the  interface  in  that  realization  and  the  average  time  at  which 
the  leading  interface  of  the  spot  appeared  at  the  particular  location  in 
the  flow.  As  a  result  of  the  alignment  process,  the  rate  of  increase  in 
veloeity  after  the  arrival  of  the  leading  interface  is  much  stronger  than 
the  rate  estimated  from  the  simple  ensemble-averaged  data,  while  in  the  in¬ 
terior  of  the  spot  the  aligned  and  the  non-aligned  ensemble  averages  are 
identical.  The  top  trace  in  figure  D.5.2  represents  the  absolute  value  of 
the  difference  between  the  two  methods  of  averaging.  It  is  obvious  that  in 
the  vicinity  of  the  leading  interface  the  difference  is  significant.  The 
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second  trace  from  the  top  represents  the  conventional  r.e.m.s  distribution 
of  the  turbulent  fluctuations  near  the  leading  interface  of  the  spot.  The 
apparent  increase  in  the  intensity  of  the  fluctuations  relative  to  the  in¬ 
terior  of  the  spot  is  of  the  same  order  of  magnitude  as  the  difference 
between  the  two  ensemble  averages.  One  may  thus  conclude  that:  the  appar¬ 
ent  turbulent  activity  near  the  leading  interface  of  the  spot  is  an  arti¬ 
fact  of  the  definition  of  the  r.e.m.s  fluctuations  and  the  method  of  aver¬ 
aging.  Making  the  same  comparison  near  the  trailing  interface  of  the  spot 
leads  to  a  similar  conclusion.  These  results  do  not  imply  that  there  is  no 
enhanced  turbulent  activity  near  the  interfaces  of  the  spot,  but  rather 
that  the  conventional  data  acquisition  and  the  definition  of  intensity  of 
the  fluctuations  is  inadequate  in  this  region.  The  time  scale  of  the 
•enhanced’  turbulent  activity  near  the  leading  interface  of  the  spot  cor¬ 
responds  to  the  Jitter  in  the  time  of  arrival  of  the  leading  interface. 
This  was  determined  from  the  histograms  of  the  arrival  times  shown  din  fig¬ 
ure  D.5.3a.  Near  the  surface,  for  example,  the  average  time  of  arrival  of 
the  leading  interface  from  the  initiation  of  counting  was  equal  to  19  msec. 
The  standard  deviation  from  the  mean  is  1.5  msec.  This  number  (also  plot¬ 
ted  in  figure  D.5.2)  indicates  that  the  duration  of  the  apparently  enhanced 
activity  corresponds  indeed  to  the  standard  deviation  in  the  time  of  arri¬ 
val  of  the  leading  interface  of  the  spot  at  the  chosen  location.  Although 
the  Jitter  in  this  time  increases  at  large  distances  from  the  surface,  the 
apparent  error  in  the  determination  of  the  r.e.m.s  values  of  the  velocity 
fluctuations  may  not  necessarily  increase,  because  only  near  the  surface 
the  arrival  of  the  spot  is  marked  by  a  significant  increase  in  velocity. 


The  duration  of  the  spot  ean  be  computed  for  each  realization  by  know¬ 
ing  the  time  of  arrival  of  the  leading  and  trailing  interface  respectively. 
A  histogram  showing  the  duration  of  the  spot  is  plotted  in  figure  D.5.3b. 
The  mean  duration  of  the  spot  is  shown  by  a  vertical  line  drawn  at  each 
level.  The  standard  deviation  (which  is  marked  by  triangles  on  the  base¬ 
line  of  each  histogram)  increases  with  increasing  distance  from  the  sur¬ 
face.  This  is  particularly  true  if  one  normalizes  the  deviation  by  the 
local  duration  of  the  spot.  The  values  obtained  are  also  shown  in  figure 
0.5.3b.  Further  examination  of  the  histograms  indicates  that  the  probabil¬ 
ity  density  distribution  is  not  centered  arround  the  mean  duration  of  the 
spot,  but  appears  to  be  lumped  in  groups  which  are  merely  represented  by 
the  mean  duration  indicated  on  the  figure.  Although  the  number  of  samples 
is  too  small  for  a  conclusive  interpretation  of  the  histograms,  the  charac¬ 
teristic  times  (shown  on  the  figure)  suggest  something  about  the  interior 
structure  of  the  spot.  If  the  spot  contains  a  finite  number  of  large  co¬ 
herent  eddies,  then  the  length  of  the  spot  depends  on  the  number  of  eddies 
within  it.  The  probability  density  distribution  for  y/ 6fc<  0.3  suggests 
that  most  of  the  spots  at  this  location  and  Reynolds  number  vary  from  one 
another  by  two  to  four  large  coherent  eddies.  The  duration  of  each  eddy  in 
this  ease  is  approximately  3  msec  and  corresponds  also  to  the  variation  in 
the  time  of  arrival  of  the  spot.  During  the  time  elapsed  between  the  arri¬ 
val  of  the  leading  interface  and  the  trailing  interface  at  the  measuring 
station  the  spot  almost  doubles  in  length.  The  variation  in  the  time  of 
arrival  of  the  trailing  interface  (fig.  D.5*3c)  is  twice  as  large  as  the 
variation  in  the  time  of  arrival  of  the  leading  interface.  Once  more  the 
length  of  the  spot  seems  to  depend  on  the  number  of  large  coherent  struc- 


Page  42 


tures  contained  In  It.  Each  structure  (Kovasznay,  Komoda  and  Vasudeva, 
1962),  Host  probably  originates  at  the  edge  of  the  laminar  boundary  layer, 
or  at  least  within  the  laminar  boundary  layer,  and  slowly  grows  in  all 
three  directions  while  propagating  downstream.  Thus  the  maximum  height  of 
the  spot  is  comparable  to  the  thiokness  of  the  local  turbulent  boundary 
layer.  Sinoe  some  structures  originate  later  than  others  and  they  also  do 
not  extend  as  far  out  as  the  edge  of  the  corresponding  turbulent  boundary 
layer,  producing  large  variations  in  the  maximum  height  of  the  spot  and  its 
length  at  the  outer  edge,  the  standard  deviation  in  the  length  of  the  spot 
at  its  outer  edge  must  also  increase  significantly  (figure  D.5.3b). 


D.6  The  Zone  Averaged  Intensities  &'  gad.  li  and.  Iha.  Zone  Average 
f  rodaftt  ji'jl'. 


In  order  to  compare  the  average  turbulent  intensity  in  an  individual 
spot  to  the  turbulent  intensity  prevailing  in  interacting  spots  or  in  a 
turbulent  boundary  layer,  one  has  to  average  the  velocity  fluctuations  dur¬ 
ing  the  time  corresponding  to  the  duration  of  the  turbulent  regions  only. 
The  method  used  presently  is  the  zone-averaging  technique  suggested  by  Ko¬ 
vasznay,  Kibens  and  Blaokwelder  (1970).  A  telegraph  signal  I(t),  corres¬ 
ponding  to  the  ooourrance  of  turbulence  and  having  an  average  value  of  : 


(17)  Y  *I*1/T)  Idt 


whioh  is  oommonly  referred  to  as  intermittenoy,  defines  a  'turbulent  zone 
average'  of  any  property  S  by  the  following: 
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T 

(18)  Ss(1/yT)  jlxSdt  . 

0 

Thus  the  turbulent  intensities  of  u  and  v  were  calculated  during  the  pas¬ 
sage  of  the  spot  as: 

T  T 

(19)  u'sjfl/yT)  [lx(U-<U»2dt]1/2  ;  v»  =  [(1/yT)  flx(V-<V»2dt]1/2 

b  0 

where,  <U>  and  <V>  are  the  ensemble-averaged  values  for  the  two  velocity 
components  U  and  V,  respectively. 

The  intensity  of  the  u  component  in  isolated  spots  is  shown  in  figure 
0.6.1,  for  the  three  measuring  stations  at  distances  of  84.5,  103  and  122 
cm  from  the  perturbation ,  together  with  the  ‘turbulent  zone  average'  in  a 
fully  turbulent  boundary  layer.  The  intensity  u’  (ordinate)  is  normalized 
by  the  free  stream  velocity,  while  distance  from  the  surface  (abscissa)  is 
normalized  by  the  height  of  a  fUlly  turbulent  boundary  layer  existing 
at  the  particular  location  and  Reynolds  number.  The  zone-averaged  u'  in 
the  spot  is  higher  by  approximately  30f  than  the  corresponding  value  in  a 
fully  developed  turbulent  boundary  layer.  The  distribution  of  u'  in  both 
flows  is  similar  for  y/6^<0.6.  It  should  be  remembered  that  the  turbulent 

w 

bulges  in  a  boundary  layer  extend  beyond  6fc,  therefore,  if  the  maximum  he¬ 
ight  of  the  spot  and  a  turbulent  bulge  are  to  be  compared,  the  spot  data 
should  be  normalized  by  its  maximum  height,  while  the  boundary  layer  data 
should  be  normalized  by  1.2  (i.e.  the  maximum  extent  of  the  turbulent 

bulges).  It  turns  out  that  this  normalization  will  not  improve  the  agree¬ 
ment  between  the  data  taken  in  isolated  spots  and  in  the  turbulent  boundary 
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layer,  because  the  maximum  height  of  the  spot  Is  also  approximately  1 . 2  4 . 
The  present  results  in  the  turbulent  boundary  layer  are  compared  with  the 
data  of  Kovasznay  et  al.  (1970).  The  agreement  between  the  two  sets  of 
data  (fig.  D.6.1)  is  very  good  in  spite  of  the  different  Reynolds  numbers 
(Refi^)  at  which  the  measurements  were  made  (13,200  and  27,500,  respective¬ 
ly).  Similar  data  for  the  normal  oomponent  of  the  fluctuations  is  shown  in 
figure  0.6.2;  It  indicates  that  the  v*  intensity  in  turbulent  spots  is 
lower  than  in  the  turbulent  boundary  layer.  Thus,  while  the  intensity  u* 
within  the  spot  is  larger  than  the  intensity  prevailing  in  a  turbulent 
boundary  layer,  the  intensity  v'  within  the  spot  is  less  than  in  the  turbu- 
lent  boundary  layer.  The  zone-averaged  -u*  v'  within  the  spot  is  also  less 
than  in  the  corresponding  fully  turbulent  boundary  layer  (figure  D.6.3.). 

0.7  Ihs.  Distribution,  al  JJ&&  tomalfla  atesMaa  an4  XaElMtonfc 

intaaaiUaa  Along  10a  Scat* 

The  zone-averaged  data  oontains  no  information  about  the  distribution 
of  physical  properties  along  the  spot.  It  would  be  naive  to  assume  that 
the  distribution  of  the  Reynolds  stress  and  turbulent  intensity  along  the 
spot  is  homogeneous,  thus  instantaneous  and  ensemble-averaged  results  are 
examined.  The  generation  of  turbulent  energy  is  given,  on  the  average,  by 
the  product  t u* v ' ( 3U/ 3y ) ]  where  u'v'  and  90/3y  refer  to  time-mean  data.  It 
is  not  easy  to  correlate  the  time  averaged  production  to  the  local  instan¬ 
taneous  occurrence  in  a  turbulent  boundary  layer.  The  instantaneous  velo¬ 
city  gradient  3U/3y  may  not  necessarily  resemble  the  mean  gradient. 
Furthermore,  the  instantaneous  velocity  gradients  in  other  directions  may 
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be  significant  in  spite  of  the  fact  that  they  vanish  in  the  mean;  thus, 
when  coupled  with  the  appropriate  stress  they  too  may  cause  a  local  produc¬ 
tion  of  turbulence.  Kline  (1978)  reviews  some  of  the  available  results  on 
the  structure  of  the  turbulent  boundary  layer  obtained  during  the  past  de¬ 
cade  and  suggests  that  70$  of  the  turbulent  production  occurs  during  the 
outward  ejection  of  low  momentum  fluid  from  the  wall  region  (bursting). 
The  rest  of  the  positive  contribution  to  production  of  turbulence  is  relat¬ 
ed  to  the  inflow  of  high  momentum  fluid  towards  the  wall  layers  (sweep). 
The  scales  of  both  phenomena  are  related  to  the  large  coherent  eddies 
(Laufer  and  Narayanan,  1970).  Since  the  mean  production  term  is  a  product 
of  two  averages,  u’v*  and  9U/ 9y,  one  may  consider  one  aspect  contributing 
to  the  production  by  assuming  9U/9y  as  being  constant  with  time  and  looking 
for  the  intermittency  in  uv;  where  uv  is  defined  by: 

(20)  uv=(U-<U»x(V-<V» 

<U>  and  <V>  are  the  same  as  in  (19). 

Typical  time  histories  of  the  instantaneous  uv  product  at  given  y  locations 
are  shown  in  figures  D.7.1a  and  b  for  two  different  realizations.  At 
y=0.12  6  most  of  the  contribution  to  uv  occurs  intermittently  along  the 

w 

spot  and  the  amplitude  of  the  instantaneous  uv  excursions  exceed  by  a  fac¬ 
tor  of  fifteen  the  long  time-average  observed  in  a  turbulent  boundary 
layer.  At  y/  6  0.21,  0.31  and  0.41  the  instantaneous  excursions  in  uv  are 

w 

often  smaller. 


A  probability  density  of  uv/u’v’  >  15,  obtained  from  90  realizations, 
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7^7 

is  shown  on  the  lower  part  of  figure  D.7.2.  (u* v' is  the  zone  averaged  uv 
product  measured  in  a  turbulent  boundary  layer  at  the  same  location  and 
same  ambient  velocity).  The  highest  concentration  of  the  spikes  exceeding 
the  prescribed  threshold  occurs  in  the  wall  region  near  the  leading  edge  of 
the  spot.  One  may  suspect  that  the  activity  near  the  interface  results 
from  the  statistical  method  used  in  obtaining  ensemble  averaged  <uv>,  (see 
the  discussion  in  section  D.5)  but  the  time  scale  of  the  enhanced  <uv>  ac¬ 
tivity  is  much  too  long  to  be  attributed  to  differences  in  the  time  of  ar¬ 
rival  of  the  spot.  The  duration  of  the  standard  deviation  in  the  spot  ar¬ 
rival  times,  which  is  also  marked  on  the  figure,  is  only  1.5  msec  while  the 
enhanced<uv>  activity  near  the  wall  lasts  10  msec. 

The  second  term  contributing  to  the  local  production  of  turbulence  is 
the  instantaneous  gradient,  9U/3 y ,  along  the  spot.  Unfortunately  the 
number  of  Instantaneously  measured  velocities  in  a  profile  does  not  suffice 
for  proper  differentiation  and  therefore  the  ensemble-averaged  data  is  pre¬ 
sented  instead  (figure  D.7.2,  top).  The  dark  lines  represent  a  low-pass 
filtered  velocity  gradient  3<U>/3y  which  was  averaged  over  100  events.  The 
dotted  lines  represent  a  portion  of  the  velocity  profiles  shown  in  figure 
D.4.1;  the  boundaries  of  the  spot  (i.e  the  2%  U^  perturbation  contour)  are 
shown  on  this  figure  for  ease  of  comparison.  The  distribution  of  the  velo- 
oity  gradient  before  the  arrival  of  the  spot  is  identical  to  the  gradient 
prevailing  in  a  laminar  boundary  layer.  Soon  after  the  arrival  of  the  spot 
there  is  an  abrupt  increase  of  approximately  40$  in  the  maximum 
value, 3<U>/3y  occurring  near  the  surface.  The  distribution  of  3<U>/3y  re- 
aaina  similar  during  the  passage  of  the  spot,  but  the  maximum  value  of 
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S<U>/aj  at  the  trailing  edge  exceeds  by  a  factor  of  2  the  maximum  value  in 
the  laminar  boundary  layer.  The  data  was  obtained  by  a  rake  of  wires  fixed 
in  space,  thus  when  3<U>/3y  near  the  wall  is  large,  there  are  too  few  data 
points  to  determine  its  value  accurately.  The  distribution  of  3<U>/3y  in 
the  central  region  of  the  spot  (marked  on  the  figure)  is  almost  identical 
to  the  distribution  of  3<U>/3y  in  a  fully  turbulent  boundary  layer. 

One  may  examine  the  contribution  to  the  production  of  turbulence  at 
each  section  along  the  spot,  by  calculating  the  value  of  <uv>x3<U>/3y.  The 
closest  location  to  the  surface  at  which  u'v'  data  is  available  is:  y=0.12 
at  this  looation  ensemble  averaged  velocity  gradient  during  the  passage 
of  the  spot  is  increased  by  a  faotor  of  1.6.  The  ratio  between  <uv>  at  the 
leading  edge  and  <uv>  at  the  trailing  edge  of  the  spot  was  found  to  be  2.5, 
therefore  the  production  at  the  leading  edge  exceeds  the  production  at  the 
trailing  edge  by  a  factor  of  1.56  (whioh  is  2. 5/1. 6).  Closer  to  the  sur¬ 
face,  however,  the  distribution  of  the  ensemble  averaged  turbulence  produc¬ 
tion  along  the  spot  may  be  different. 


CHAPTER  E 


b.  i .  Zb&  Hrib  Xalaauz  Cisll 


In  spit*  of  their  complexity,  isolated  spots  may  still  be  regarded  as 
fundMiental  struotures  in  the  transition  process.  In  this  experiment  they 
were  generated  artificially  by  an  eleotrio  discharge  or  a  Jet  emanating 
momentarily  from  a  small  orifice  in  the  plate.  An  insertion  of  a  roughness 
element  instead  of  a  momentary  disturb? 'ice  results  in  a  generation  of  a 


continuous  cone  of  turbulence  spreading  laterally  with  increasing  distance 
downstream.  The  cone  of  turbulence  degenerates  into  an  array  of  spots  whe¬ 


never  the  Reynolds  number  is  reduced  below  a  certain  threshold  level.  A 
further  reduction  in  Re  does  not  trip  the  boundary  layer.  In  the  transi¬ 
tional  Reynolds  number  range,  the  spots  appear  randomly  in  time  although 
they  oeour  directly  downstream  of  the  perturbance.  In  order  to  simplify 
the  data  aquisition  process  and  the  ooneoamitant  statistical  analysis,  one 
may  produce  spots  at  regular  intervals,  from  a  single  source  and  let  them 
interact  downstream.  The  interaction  process  generates  a  quasi  steady  tur¬ 
bulent  boundary  layer  whioh  is  disoussed  in  this  chapter.  The  spots  are 
considered  merged  whenever  the  intermittenoy  faotor  becomes  unity  for  y/H  < 
0.3  (sinoe  H  corresponds  to  1.2  in  a  fully  turbulent  boundary  layer  in 
absenoe  of  pressure  gradient  and  the  latter  is  Intermittent  whenever  y/6fc  > 
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The  present  experiment  could  be  perforated  most  naturally  by  generating 
a  continuous  array  of  spots  at  equal  time  intervals.  In  such  an  experiment 
the  data  would  be  acquired  continuously ,  provided  one  channel  was  connected 
to  the  spot-generator  for  time  reference.  However,  the  spark  perturbing 
the  laminar  boundary  layer  is  generated  by  a  high  voltage  source  (20  KV) 
which  induces  spurious  signals  on  the  data  acquisition  lines.  It  is  thus 
advantageous  to  digitize  the  data  while  the  spark  generator  is  idle.  By 
interrupting  the  electric  discharge  during  the  acquisition  of  data,  one 
generates  a  train  of  spots  whose  length  depends  on:  the  frequency  of  the 
disturbance,  the  distance  between  the  generator  and  the  measuring  station 
and  U  .  Only  the  first  and  the  last  spots  in  the  train  are  different  from 

CO 

the  rest  because  they  interact,  at  one  of  their  interfaces,  with  the  unper¬ 
turbed  laminar  boundary  layer  rather  than  a  neighbouring  spot;  this  is  ad¬ 
vantageous  because  it  enables  one  to  assess  the  effects  of  partial  interac¬ 
tion  on  the  overall  structure  of  the  spot. 

An  array  of  five  successive  spots  originating  from  a  point  source  at 
predetermined  time  intervals  was  monitored  by  a  rake  of  hot-wires  located 
on  the  plane  of  symmetry  at  several  distances  downstream  of  the  distur¬ 
bance.  The  degree  of  interaction  between  adjacent  spots  could  be  altered 
by  either  changing  the  frequency  of  their  generation,  or  by  changing  the 
downstream  distance  at  which  the  spots  were  monitored.  Two  sets  of  veloci¬ 
ty  histories  are  shown  in  figures  B. 1.1a  and  b.  Individual  spots  are 
dearly  recognizable  in  figure  E.l.la,  where  the  time  interval  between  ad¬ 
jacent  spots  was  66  milliseconds  and  the  measurement  was  made  86  cm  down¬ 
stream  of  the  disturbance;  however,  when  the  time  interval  was  decreased 
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to  >10  Billiseconds  at  Xs»124  as,  a  continuously  turbulent  velocity  signa¬ 
ture  was  observed  over  a  large  fraotion  of  the  boundary  layer  thickness. 
This  signal  is  quite  similar  to  velocity  histories  aeasured  in  a  fully  tur¬ 
bulent  boundary  layer  although  a  periodic  behaviour  is  visible  at  large 
distances  from  the  surface. 

A  velocity  profile  obtained  by  ensemble-averaging  the  data  over  all 

the  events  and  then  averaging  over  time  (excluding  the  first  spot  and  the 

ealm  region  behind  the  last  spot  in  an  array)  is  shown  in  figure  E.1.2  for 

the  various  oases  considered.  Measurements  made  in  the  fully  turbulent 

boundary  layer  at  the  same  location  on  the  plate  are  marked  by  diamonds  on 

the  figure.  The  interaction  of  suooesalve  spots  produces  a  mean  velocity 

profile  which  tends  to  a  conventional  profile  as  the  interaction  becomes 

stronger.  An  array  of  transitional  spots  will  give  rise  to  a  logarithmic 

profile  whenever  T  =  ATU  /Xs  <  0.4  ,  where  AT  is  the  time  interval  between 
n  <» 

adjacent  spots  ( Wygnanski  1978). 

When  the  velooity  histories  are  ensemble  averaged  conditionally  to  the 
first  spark,  the  periodio  behaviour  introduced  by  the  disturbance  becomes 
apparent.  Contours  of  constant  velooity  patterns  are  plotted  in  figures 
E.1.3a,b  and  o.  The  oontours  representing  excess  of  velocity  relative  to 
the  mean  are  shown  as  solid  lines,  while  the  oontours  showing  defect  are 
marked  by  dashed  lines.  The  perturbation  contours  in  figures  B.1.3  refer 
to  the  laminar  velooity  profile  existing  at  the  measuring  station  in  ab¬ 
sence  of  any  disturbance.  When  the  dimensionless  interval  between  adjacent 
spots  is  T  =1.16,  the  velooity  perturbation  contours  (fig.  E.1.3a)  are 
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very  similar  to  those  plotted  for  the  isolated  spots  (fig.  D.1.1).  The 
contours  showing  positive  velcoity  perturbation  which  stretch  over  a  long 
time  period  behind  an  isolated  spot  appear  to  join  together  in  figure 
B. 1.3b,  where  Tn=0.465.  This  time  interval  is  sufficiently  short  to  result 
in  a  continuously  turbulent  signal  near  the  surface.  Two  velocity  profiles 
are  drawn  to  the  same  scale  in  the  upper  left  corner  of  figure  E.I.30,  one 
shows  a  velocity  distribution  in  laminar  flow,  the  other  the  turbulent  vel¬ 
ocity  profile,  both  existing  at  Xss124  cm.  The  perturbation  contours  shown 
in  the  figure  correspond  approximately  to  the  difference  between  the  two 
profiles  when  the  Blasius  profile  is  taken  as  reference.  The  sequential 
interaction  of  transitionsal  spots  contains  therefore  some  of  the  most  im¬ 
portant  elements  of  the  turbulent  boundary  layer  which  Justifies  this  some¬ 
what  artificial  experiment. 

E.2  i&a  £aagafel9=A3taEa&sA  acaatocft  sol  JgRtealnmsnt 
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In  this  section  the  ensemble  averaged  structure  resulting  from  the  in¬ 
teraction  of  spots  originating  from  a  single  source  is  examined.  For  this 
purpose  the  shortest  time  interval  (27  msec)  between  successive  spots  was 

chosen  in  order  to  assure  that  the  dimensionless  time  (T  )  at  the  first 

n 

measuring  station  (Xss84.5  cm)  was  less  then  0.4.  As  mentioned  before, 
(sec.  E. 1),  Tn*0.4  represents  a  threshold  value  below  which  the  spots 
merge  to  form  a  continuous  region  of  turbulent  flow.  Thus,  for  the  chosen 
time  interval,  the  free  stream  velocity  and  distances  from  the  source,  the 
value  of  Tn  varied  from  0.31  to  0.22.  The  mean  velocity  perturbation  con¬ 
tours  on  the  plane  of  symmetry,  103  cm  downstream  from  the  disturbance  are 
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represented  in  figure  E.2.1.  The  perturbation  velocity  was  calculated  re¬ 
lative  to  the  laminar  boundary  layer  profile  (using  eq.  9),  repeating  the 
procedure  adopted  for  the  isolated  spots  (fig.  D. 1.1).  The  first  spot  in 
the  array  maintains  a  maximum  height  of  H/ 6^=1. 2,  which  is  equal  to  the 
maximum  height  attained  by  isolated  spots  in  spite  of  the  fact  that  its 
trailing  interface  was  overtaken  by  the  leading  interface  of  the  following 
spot.  All  spots,  with  the  exception  of  the  last  one,  lose  the  ’calm’  high 
velocity  region  which  always  follows  an  isolated  spot.  As  a  result  of  the 
interaction,  the  absolute  level  of  the  perturbation  (U-U^J/U^  is  reduced  by 
approximately  20%  in  comparison  with  similar  perturbations  observed  in  the 
isolated  spots. 

The  —2J6  perturbation  contour,  assumed  to  represent  the  average  posi¬ 
tion  of  the  interface,  oscillates  around  a  mean  location  y/6t=0.78.  The 
convection  velocity  of  the  leading  interface  is  fairly  constant  (0.88Ura  ) 
and  is  independent  of  the  degree  of  interaction  between  successive  spots 
which  is  determined  by  "Tn".  This  is  in  agreement  with  the  data  of  Savas 
(1979)  who  observed  that  the  celerity  of  interacting  spots  is  constant 
(0.881^  ),  independent  of  their  scale.  The  celerity  of  the  trailing  inter¬ 
face,  however,  increases  from  0.55U  at  T  =0.96  to  a  value  of  about 

00  n 

0.80^  at  Tn=  0.26.  With  the  exception  of  the  extreme  spots  in  the  train, 
the  perturbation  contours  of  the  streamwise  component  of  velocity  can  be 
divided  into  two  families.  In  the  outer  region  of  the  boundary  layer  the 
defect  contours  are  wavy  and  parallel  to  one  another.  Near  the  surface  the 
excess- velocity  contours  are  almost  straight  and  parallel  to  the  solid  sur¬ 
face.  The  crests  of  the  defect  contours  correspond  to  the  high  points  on 
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the  spots  while  the  valleys  correspond  to  the  potential  fluid  engulfed  in 
the  region  of  interaction  between  adjacent  spots.  One  may  draw  straight 
lines  (fig.  E.2.1)  joining  the  bottom  of  the  valleys  or  the  crests  of  all 
contours  for  a  given  ensemble  averaged  case.  The  inclination  of  these 
lines  to  the  surface  is  15°  for  the  valleys  and  10°  for  the  orests  when 
based  cxi  an  average  convection  velocity  0.88(1^.  It  is  interesting  to  com¬ 
pare  these  results  with  the  data  presented  by  Brown  and  Thomas  (1977)  who 
correlated  the  shear-stress  near  the  surface  with  the  streamwise  velocity 
at  four  locations  in  the  turbulent  boundary  layer  over  short-time  periods. 
The  location  of  the  probes  at  which  the  correlation  was  maximum  at  TsO 
oould  be  nearly  connected  by  a  straight  line  inolined  at  18°  to  the  wall. 
The  large  ooherent  eddies  in  a  fully  turbulent  boundary  layer  when  visual¬ 
ized  by  smoke  are  inclined  to  the  surface  at  approximately  the  same  angle 
(see  Falco  1977).  The  turbulent  non- turbulent  interface  for  the  interact¬ 
ing  spots  was  detected  in  the  usual  manner  (e.g  WSF  1975)  for  the  case  when 
Tn*0.465.  The  average  shape  of  the  interface  is  clearly  distinguishable 
between  y/<5j*1.3  and  y/6^2.2  (fig.  E.  1.3b).  The  leading  interface  is  in¬ 
clined  to  the  surface  of  the  plate  at  18°.  Since  the  orests  of  the  spots 
or  the  valleys  interconnecting  them  are  not  always  at  the  same  elevation 
above  the  surface,  it  is  difficult  to  define  the  average  position  of  the 
interface  for  y/6^  >2.2  and  y/6^  <3.2  unambiguously.  One  may  notice  in 
figure  E.1.3b  that  the  2%  defect  perturbation  contour  no  longer  coincides 
with  the  looation  of  the  interface.  When  Tn  <  0.4,  in  fact,  no  contour  re¬ 
presenting  oonstant  velocity  perturbation  is  even  approximately  parallel  to 


the  interface 


The  flow  field  within  the  region  of  interaction  may  be  examined  by 
plotting  the  ensemble  averaged  perturbation  velocity  vectors  in  the  X-Y 
plane  obtained  from  a  rake  of  X-wire  probes  (fig.  E.2.2).  The  perturba¬ 
tions  considered  are  relative  to  the  laminar  boundary  layer  (using  eq.  9 

and  10  for  <0  >  and  <V  >).  The  -2 i  U  perturbation  contour  is  also  replot- 
P  P  00 

ted  in  figure  E.2.2  for  reference.  A  negative  value  of  the  normal  velocity 
component  is  observed  only  once  near  the  leading  interface  of  the  first 
spot,  where  the  flow  is  free  from  interactions  with  other  turbulent  struc¬ 
tures.  The  normal  velocity  component  during  the  passage  of  the  other  spots 
in  the  array  is  rarely  negative,  resembling  the  prevailing  direction  of  the 
flow  inside  a  transitional  spot  (fig.  D.1.2a).  There  is  however  an  appar¬ 
ent  periodicity  in  the  magnitude  of  the  outflow  from  the  surface,  matching 
the  time  interval  between  spots. 

The  rate  of  entrainment  of  irrotational  fluid  by  a  turbulent  boundary 
layer  is  always  of  interest  because  it  is  coupled  with  the  growth  of  the 
boundary  layer  and  often  serves  for  calculation  purposes  (Head  1958).  The 
entrainment  of  irrotational  fluid  by  an  array  of  successive  spots  on  their 
plane  of  symmetry  may  not  only  give  us  some  insight  into  the  kinematic  des¬ 
cription  of  the  entrainment  process  but  also  into  the  relative  importance 
of  some  three  dimensional  effects.  The  following  calculations  were  done  in 
the  same  manner  as  for  the  isolated  spot  (section  D.2): 

(i)  The  entrainment  was  calculated  relative  to  the  -2% 

U  contour. 

oo 

(ii)  The  velocity  field  in  the  X-Y  plane  is  provided  by  the 


rake  of  X-wires. 


.1 
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(iii)  The  st ream wise  component  of  the  eelerity  of  the  inter¬ 
face  (U-U^s -2t  Uo)  is  obtained  from  the  arrival  times 
of  the  interfaces  at  the  three  x-stations  at  which 
measurements  were  made.  The  celerities  of  the  leading 
and  trailing  interfaces  far  away  from  the  surface  are 
0.90c  and  0.8U,,,  respectively.  It  is  difficult  to  esti¬ 
mate  the  celerity  of  the  interface  precisely,  in  par¬ 
ticular  in  the  region  where  two  adjaoent  spots  merge, 
(i.e.  at  the  bottom  of  the  turbulent  valleys). 
Fortunately  however,  the  difference  between  the  celer¬ 
ities  of  the  leading  and  trailing  interface  is  not  as 
large  as  in  the  isolated  spot,  and  the  possible  error 
inourred  in  the  calculations  of  the  entrainment  is  not 
large . 

(iv)  The  normal  component  of  the  eelerity  of  the  interface 
is  estimated  to  be  a  constant  independent  of  y: 
Vint/UO0  =6.5x10  .  This  estimate  is  based  on  the  rate 
of  growth  of  the  apex  of  a  typical  spot  in  the  array 
with  downstream  distance.  In  isolated  spots  V^nt  var¬ 
ies  with  normal  distanoe  from  the  "overhang"  and  al¬ 
though  the  same  variation  may  still  exist  in  the  pre¬ 
sent  case  the  total  variation  in  the  height  of  the  in¬ 
terface  (which  is  approximately  15%  6.)  does  not  war- 
rant  the  added  difficulties  in  calculating  the  entrain¬ 
ment. 

(v)  The  total  amount  of  fluid  entrained  by  a  typical  spot 


in  an  array  on  its  plane  of  syoaietry  is  obtained  by  in¬ 
tegrating  the  nonaal  velocity  component  relative  to  the 
interfaoe  over  the  length  of  the  interface. 

For  the  specific  oase  chosen,  the  rate  of  entrainment  is  6.5x10“^ 

2 

a  /sec;  40?  of  whioh  occurs  through  the  leading  interfaoe. 

The  flow  field  relative  to  the  ensemble  averaged  interfaoe  of  a  typi¬ 
cal  front  in  an  array,  (fig.  E.2.3a)  nay  be  compared  qualitatively  with 
the  flow  field  relative  to  a  bulge  in  turbulent  boundary  layer,  observed  by 
Blaokwelder  (1970)  and  reproduoed  in  figure  E.2.3b.  The  aoale*  of  the 
bulge  and  the  spot  are  different,  because  the  latter  oontains  a  number  of 
eddies  ( section  E.3)  whioh  do  not  retain  their  precise  phase  relative,  to 
the  interface;  ..  nevertheless  the  general  flow  field  relative  to  the  in¬ 
terfaoe  is  similar  in  both  oases. 

One  may  compare  the  entrainment  calculated,  by  assuming  that  the  spots 
in  the  array  are  two  dimensional,  with  the  entrainment  oaloulated  from 
their  overall  growth-rate;  or,  the  gross  behaviour  of  the  oomparable  two 
dimensional  turbulent  boundary  layer.  The  rate  of  entrainment  of  lrrota- 
tional  fluid  in  a  turbulent  boundary  layer  is  given  by: 

(21)  dq/dX*d/dX{U  (6.-6.*))  (Head  1958) 


and  is  tabulated  below: 
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Table  4 


Uoo(6-6*)x103[ 

m2/seo] 

O/U.JxCdq/dXjxlO4 

((1/Uoo)xq/UooAT)x10 

X9  [cm.] 

SPOT  TRAIN 

T.B.L 

SPOT  TRAIN 

T.B.L 

2-D  CALCULATIONS 

84.5 

144.7 

143.4 

139 

131 

— 

103 

170.5 

167.6 

133 

132 

27 

122 

195.7 

192.8 

... 

___ 

___ 

The  rate  of  growth  of  spots  In  a  train  is  similar  to  the  rate  of  growth  of 
the  turbulent  boimdary  layer  and  is  roughly  equal  to  1 . 3%  of  the  free  stre¬ 
am  velocity.  The  entrainment  velocity  calculated  assuming  that  the  flow  is 
2-dimensional  is  5  times  smaller I  The  discrepancy  between  the  two  methods 
of  calculation  stresses  the  relative  importance  of  the  entrainment  in  the 
X-Z  plane  (i.e  by  the  V  and  W  components  of  velocity).  Furthermore,  calcu¬ 
lations  based  on  ensemble  averaged  data  underestimate  the  entrainment  be¬ 
cause  both  the  contortions  of  the  interface  and  the  velocity  relative  to  it 
are  smoothed  by  the  averaging  process. 


e.3  A  single  RqaUaaUon 


One  often  asks  how  much  information  is  lost  by  the  averaging  process 
and  how  auoh  distortion  may  be  introduced  by  a  particular  eduction  scheme. 
A  simple  answer  may  be  obtained  by  comparing  an  individual  realization  with 
the  ensemble-averaged  data.  The  velocity  history,  measured  with  a  rake  of 
normal  wires,  was  low-pass  filtered  at  250  Hz  (fig.  E. 3- 1 )  and  used  for 
plotting  the  streamwise  perturbation  contours  (fig.  E. 3.2) .  A  comparison 
is  made  with  the  3%  ensemble  averaged  perturbation  contours  which  are  re¬ 
presented  in  the  figure  by  the  dashed  lines. 

Pour  successive  spots  oan  be  identified  by  the  outward  excursion  of 
the  -3$  contour,  whiah  ocours  at  the  same  time  as  the  corresponding  ensem¬ 
ble  averaged  results.  However,  while  the  ensemble  average  contours  show 
only  4  eddies  corresponding  to  the  four  artificially  evoked  spots,  the  sin¬ 
gle  realization  is  composed  of  at  least  twice  as  many  eddies.  One  may 
identify  2  eddies  per  spot  which  are  marked  by  A  and  B  on  the  figure.  With 
the  exception  of  the  leading  spot,  the  length  scale  of  the  average  struc¬ 
ture  is  approximately  66  which  is,  at  least,  comparable  to  the  length  aoale 
of  the  large  coherent  eddies  eduoed  by  correlation  techniques  (Kovasznay  et 
al.  1970),  or  to  the  structure  evoked  by  letting  a  spot  interact  with  a 
tripped  turbulent  boundary  layer.  The  number  of  structures  observed  in  a 
single  realization  of  an  isolated  spot  at  oomparable  Re  is  muoh  larger  than 
two.  It  is  thus  suggested  that  the  number  of  large  eddies  in  a  spot  is  re- 
duoed  by  the  interaction  process;  this  is  accompanied  by  a  corresponding 
reduction  in  the  scale  of  the  spot  and  also  a  reduction  in  its  capability 
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to  engulf  non-turbulent  flow. 

The  flow  in  a  single  realization  was  also  examined  by  plotting  the  in¬ 
stantaneous  velocity  vectors,  relative  to  the  convection  velocities  of  the 
interface  in  the  T-I  plane  (fig.  E. 3. 3) .  (i.e  0.881^  was  subtracted  from 
the  measured  streamwise  component  of  velocity,  and  an  average  normal  compo¬ 
nent  of  the  interface  celerity  (0.0065U  )  was  subtracted  from  the  normal 

oo 

component  of  velocity).  The  approximate  length  of  the  eddies,  which  appear 
in  this  figure,  is  only  26  as  was  deduoed  from  multiplying  their  duration 
time  (3-4  msec)  by  the  streamwise  convection  velocity  (0.88U  ).  Thus  it 

OO 

appears  that  the  scale  of  the  large  eddies  in  interacting  transitional 
spots  is  comparable  to  the  large  scales  in  a  normal  turbulent  boundary 
layer. 


e.4  Xha  BiatrttMUQn  a £  ffean  tolaaitx. 

The  evolution  of  the  mean  profile  in  an  array  of  spots  is  examined  in 

figures  E.4.1  and  E.4.2.  In  figure  E.4.1,  the  ensemble  averaged  streamwise 

velooity  profiles  taken  at  8  msec  intervals  (solid  lines)  are  compared  with 

the  laminar  velocity  profile  (dots),  existing  at  the  same  location  (X  =103 

s 

cm).  The  ensemble-averaged  profiles  at  each  section  are  also  compared  with 
the  time  averaged  profiles  along  three  central  spots  in  the  array  (figure 
E.4.2).  The  appropriate  perturbation  velocity  contours  are  replotted  in 
figures  E.4.1  and  E.4.2  for  ease  of  comparison.  The  ensemble-averaged  and 
the  time-averaged  profiles  are  almost  identical,  provided  the  first  spot  in 
the  trains  is  exoluded  from  the  comparison.  The  velocity  profile  in  the 
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latter  aay  deviate  by  O.IU^  from  the  time  average  profile  In  the  Interact¬ 
ing  spots. 

The  ensemble  averaged  profiles  plotted  on  a  semilogarithmic  scale  and 
taken  at  three  successive  locations  ranging  from  84.5-122  cm  downstream 
from  the  perturbation  are  plotted  in  figure  E.4.3  (upper  part).  The  Rey¬ 
nolds  number  based  on  the  momentum  thickness  increases,  with  increasing 
downstream  distance  from  1200  to  1700,  and  the  wake  component  increases 
from  1.5  to  2.7.  When  comparing  these  results  with  the  corresponding  lo¬ 
garithmic  profiles  in  the  turbulent  boundary  layer  at  the  same  Reynolds 
numbers  (figure  E.4.3,  lower  part),  one  may  observe  that: 

(i)  The  velocity  distribution  in  the  wall  region  converges 

to  the  same  logarithmio  law,  but  the  friction  velocity 
in  the  case  of  the  interacting  spots  is  4%  higher  than 
in  the  fully  developed  turbulent  boundary  layer. 

(ii)  The  wake  component  of  the  interacting  spots  is  smaller 

by  40?  at  X  =  84.5  cm  and  by  10$  at  X  =  122  on. 
s  s 

One  may  attribute  the  high  friction  velocity  and  the  low  wake  compo¬ 
nent  seen  in  interacting  spots  to  a  smaller  number  of  coherent  eddies  con¬ 
tained  in  the  spot  relative  to  the  number  of  eddies  contained  in  a  turbu¬ 
lent  boundary  layer  during  the  same  period  of  time.  Similar  behaviour  was 
observed  in  a  turbulent  boundary  layer  with  increasing  x,  i.e,  the  friction 
velooity  decreases  while  the  momentum  thickness  increases.  Seven  of  the 
ensemble  averaged  profiles  (measured  at  X  *103  cm)  marked  by  arrows  in  fig- 
ure  E.4.1  are  replotted  on  semilogarithm io  seale  (fig.  E.4.4)  and  the  con- 


stants  derived  from  fchls  plot  are  ooapared  with  the  universal  constants  re¬ 
corded  over  the  years  in  a  turbulent  boundary  layer. 

The  values  of  the  friction  velocity  and  of  the  wake  component  for 
these  profiles  are  tabulated  below: 


Table  5 

Time  [msec] 

au/ut 

U  /U 

*[  CO 

64 

2.05 

0.0435 

88 

2.27 

0.0435 

92 

2.73 

0.043 

96 

2.50 

0.043 

100 

2.10 

0.043 

104 

1.40 

0.044 

108 

1.80 

0.044 

T.B.L 

2.70 

0.0415 

The  results  indieate  that  the  miction  velocity  in  the  central  region 
of  the  spot  is  higher  than  in  the  turbulent  boundary  layer  while  the  con- 
oommitant  wake  component  is  lower.  One  may  conclude  that,  for  Re,0<  1700, 
the  interaction  among  successive  spots  is  not  complete  and  that  the  distri¬ 
bution  of  the  miction  velocity  and  of  the  wake  component  along  the  spots 


is  not  homogenuous. 
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E-5  Turbulent  iat.wai.uaa  and  Baynaida  S&coaa 

The  calculation  of  turbulent  intensities  and  Reynolds  stress  during 
the  passage  of  a  train  of  spots  raises  similar  difficulties  to  the  ones  en¬ 
countered  in  the  analysis  of  isolated  spots.  The  flow  is  non- stationary, 
thus  the  analysis  of  the  data  should  be  conditioned  to  the  leading  inter¬ 
face  of  the  second  spot  and  the  duration  of  the  sampling  window  should  con¬ 
tain  an  integral  number  of  events  (ATxN  when  AT  is  the  time  interval 
between  successive  spots) .  The  fluctuations  were  calculated  by  subtracting 
the  ensemble  averaged  data  from  the  instantaneous  values. 

The  level  of  the  u  fluctuations  in  the  interacting  train  of  spots,  far 
from  the  surface,  is  lower  than  the  level  encountered  in  the  isolated  en¬ 
semble  averaged  spot  by  approximately  20%  (fig.  E.5.1).  The  regularity  at 
which  these  spots  are  generated,  which  smoothes  the  contortions  of  the 
outer  interface,  may  ad  so  cause  smaller  variability  in  the  instantaneous 
velooity  and  thus  reduce  the  apparent  intensity  of  the  turbulent  fluctua¬ 
tion.  The  intensity  in  the  vicinity  of  the  leading  interface  of  the  first 
spot  in  the  train  is  also  higher  than  the  intensity  prevailing  in  other 
spots  and  resembles  the  effect  discussed  in  section  D.5  (fig.  E.5.1).  In 
order  to  oheok  whether  this  effect  stems  from  the  jitter  in  the  time  of  ar¬ 
rival  of  the  leading  spot  at  the  measuring  station,  a  histogram  of  the  ar¬ 
rival  times  was  calculated  together  with  the  standard  deviation  from  the 
mean  time  of  arrival  (see  insert  on  fig.  E.5.2).  This  data  was  compared 
and  found  to  be  identical  to  the  standard  deviation  calculated  previously 
for  isolated  spots,  suggesting  that  the  increase  in  the  intensity  u'  near 
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the  leading  interface  is  probably  an  artifact  of  the  statistical  procedure. 
The  zone-averaged  distributions  u'  and  v*  are  shown  in  figure  E.5.3  and 
compared  to  the  zone-averaged  data  in  a  fully  turbulent  boundary  layer. 
The  agreement  between  the  two  distributions  is  good.  The  values  of  the  u' 
intensity  drops  below  the  values  observed  in  the  turbulent  boundary  at  y  > 
6  as  a  result  of  the  decrease  in  the  height  of  the  spots  within  the  array. 

The  intensity  of  the  normal  fluctuations  v'  is  somewhat  lower  for  y  > 
0.5  6  .  The  agreement  between  the  distribution  of  uv  in  interacting  spots 
and  the  turbulent  boundary  layer  is  poorer  than  between  the  intensities  of 
u  and  v  separately.  Near  the  surface  (y  <  0.4  6  )  the  value  of  -u’v'  in 

V 

the  ease  of  the  interacting  spots  is  higher  than  in  the  turbulent  boundary 
layer,  while  in  the  intermittent  region  (y  >  0.6  6fc)  it  drops  below  the 
prevailing  value  measured  in  a  turbulent  boundary  layer  by  as  much  as  3051 
locally  (fig  E.5.4).  An  examination  of  the  uv  signal  reveals  that  the 
largest  excursions  in  uv  occur  in  the  wall  region  near  the  leading  inter¬ 
face  of  each  spot  (marked  by  "LE"  in  figures  E.5.5a  and  b)  suggesting  that 
turbulent  production  does  not  occur  randomly  in  the  spot. 


5 


e.  6  Observations  Madft  Eac.  Asa i  Ecaa  laa  Elans.  &£  amatoa 


Thus  far,  the  flow  field  resulting  from  the  interaction  of  successive 
spots  generated  in  tandem,  was  examined  only  on  the  plane  of  symmetry,  di¬ 
rectly  downstream  of  the  perturbation.  Some  effects  of  the  interaction  on 
the  spanwise  evolution  of  the  spots  are  discussed  in  this  section. 

Measurements  were  made  with  a  rake  of  normal  wires  located  at  z/X  =  0.15, 
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while  a  single  wire  looated  on  the  plane  of  symmetry  (i.e  at  z/X  =  0)  moni- 

5 

tored  tiie  passage  of  the  spots. 

A  typical  velocity  history,  resulting  from  a  passage  of  an  isolated 
spot  is  shown  in  figure  E.6.1.  The  uppermost  trace  in  this  figure  was  re¬ 
corded  by  the  monitoring  wire  at  z/X  =  0  while  the  other  9  traces  refer  to 

3 

the  rate  of  wires  located  at  z/X  =  0.15  and  0.07  <  y/<5.  <  0.995  as  marked. 

s  t 

The  wave  packet  following  the  spot  is  clearly  seen  at  y/6  <  0.12  (see  WHK, 

w 

1979).  The  duration  of  the  spot  at  z/X  =0.15  is  much  shorter  than  on  the 

9 

plane  of  symmetry  as  may  be  deduced  by  comparing  the  duration  of  the  turbu¬ 
lent  signatures  at  z/X  =0  with  signatures  recorded  at  z/X  *0.15. 

8  S 

The  velooity  history  during  the  passage  of  a  train  containing  5  spots 

is  shown  in  figures  E.6.2a,b.  Although  the  trace  recorded  at  z=0  shows 

clearly  the  passage  of  5  suooessive  events  at  intervals  of  33  msec  (i.e  T 

n 

*  0.32),  the  trace  recorded  at  z/X  =  0.15  shows  only  either  3  or  4  distino- 

3 

tive  turbulent  regions  (see  fig.  E.6.2a  and  fig.  E.6.2b  respectively). 
The  second  and  fourth  spots  are  missing  in  figure  E.6.2a  while  the  third 
spot  is  missing  from  figure  E.6.2b. 

A  telegraph  signal  representing  intermittency  is  plotted  in  figure 
E.6.3  and  a  probability  density  distribution  for  a  region  to  be  turbulent 
is  plotted  in  figure  E.6.4.  The  first  spot  in  the  train  is  on  the  average 
the  widest.  The  second  spot  in  the  train  is  narrower  than  the  rest  because 
in  50*  of  the  events  the  turbulent  region  did  not  reach  z/X  =  0.15.  The 

9 

probability  distribution  indicates  that  the  remaining  3  spots  in  the  train 
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are  wider  than  z/X  =0.15  in  approximately  75$  of  the  events.  This  effect 
s 

may  be  linked  with  the  "eddy  transposition'1  phenomenon  observed  by  Savas 
(1979),  however  this  could  not  be  proven  in  the  experimental  setup  des¬ 
cribed.  The  prolonged  interaction  of  spots  in  the  streamwise  direction  may 
simply  cause  a  reduction  in  their  width. 

The  wave  packet  trailing  every  individual  spot  all  but  disappears  as  a 


result  of  the  interaction  process. 
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CHAPTER  F 

LATERAL  INimCIIQM.  fiESMEM  2HQ  SEQIS. 

The  evolution  of  two  spots,  evoked  simultaneously,  was  investigated  by 
Elder  (I960).  A  plan  view  of  the  intermittency  contours  in  the  region  of 
interaction  between  the  two  spots  indicated  that  the  spreading  rate  of  each 
spot  was  not  influenced  by  the  presence  of  its  neighbour.  In  the  first 
part  of  this  experiment  the  vertical  evolution  of  two  spots  was  examined  by 
monitoring  the  velocity  field  in  the  plane  of  symmetry  between  two  spots 
generated  simultaneously  from  two  sources  located  15  cm  apart  in  the  span- 
wise  direction.  The  ensemble  average  contours  of  the  streamwise  component 
of  the  velocity  perturbation  on  the  plane  of  symmetry  between  the  spots  is 
shown  in  figure  F. 1  (by  solid  lines). 

Switching-off  one  of  the  spots  generators  and  recording  the  data  at 

the  same  location  (i.e  z/X  =  0.06  for  a  single  spot)  provides  direct  esti- 

s 

mate  of  the  importaroe  of  such  interaction  (fig.  F.1  by  dashed  lines). 
The  effects  of  the  interaction  seem  negligible  near  the  surface  although 
the  maximum  height  of  the  contours,  produced  by  the  interacting  structures, 
may  be  somewhat  increased.  Elder's  observation  related  to  the  superposi¬ 
tion  of  spots  is  thus  valid.  The  ensemble  averaged  velocity  vectors  (figs. 
F.2a,b),  obtained  from  a  rake  of  X-wires,  indicates  that  two  large  eddies 
may  be  required  to  represent  the  velocity  field  in  the  laterally  interact¬ 
ing  spots  (i.e  the  average  sense)  while  only  one  may  suffice  to  represent 
the  velooity  field  in  isolated  spots. 


By  delaying  the  generation  of  one  spot  relative  to  the  other,  one  may 
observe  the  interaction  of  two  spots  in  a  skewed  configuration  (fig.  F.3). 
All  the  preliminary  data  available  confirms  the  superposition  observation 
of  Elder.  The  relatively  small  effects  produced  by  the  interaction  of 
neighbouring  spots  are  in  accord  with  the  supposition  that  the  spot  con¬ 
tains  a  number  of  large  eddies  which  are  continuously  interacting  among 
themselves. 

As  was  shown  before,  the  pattern  of  spots  created  in  tandem  is  not  the 
best  for  generating  a  synthetic  turbulent  boundary  layer  (in  agreement  with 
the  conclusions  of  Savas  1979).  However,  since  the  growth  of  the  spots  is 
non-sensitive  to  the  presence  of  other  spots  at  different  spanwise  loca¬ 
tions,  an  optimistic  way  of  generating  the  turbulent  boundary  layer,  plus 
marking  large  structures  in  it  would  be:  to  spread  over  the  surface  a 
small  number  of  spot  generators  which  may  be  activated  randomly. 


CHAPTER  G 


MHcmpma  amfiKa 

The  simulation  of  a  fully  developed  turbulent  boundary  layer  by  an 
array  of  artificially  generated  transitional  spots,  which  are  allowed  to 
interact  under  controlled  conditions,  was  examined.  The  turbulent  boundary 
layer  is  usually  described  statistically  using  time  averaged  quantities  of 
mean  velocity,  turbulent  intensity,  Reynolds  stress  and 
spaee-time-correlations.  More  recently  conditional  sampling  techniques  are 
used  in  an  attempt  to  understand  the  mechanism  governing  the  structure  of 
the  turbulent  boundary  layer. 

An  array  of  successive  transitional  spots,  generated  from  a  single 
source  at  regular  time  intervals,  simulates  the  mean  velocity  profile  en¬ 
countered  in  a  turbulent  boundary  layer  fairly  well.  A  comparison  of  the 
logarithmic  and  the  wake  components  of  the  velocity  profile  reveals  no 
major  differences  between  the  synthetic  and  the  natural  turbulent  boundary 
layers.  (The  term  synthetic  boundary  layer  was  coined  first  by  Coles  and 
Barker  1975). 

Although  the  periodicity  at  which  the  spots  were  generated  could  easi¬ 
ly  be  detected  by  considering  phase-locked  averaged  data,  the  turbulent  in¬ 
tensity  and  the  Reynolds  stress  (averaged  over  time)  within  the  turbulent 
zone  of  the  synthetic  and  the  natural  boundary  layers  were  very  similar. 
The  excursions  of  the  interface  in  the  synthetio  boundary  layer  were,  of 
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course,  periodic  and  limited  in  their  outward  extent  in  comparison  to  the 
natural  boundary  layer.  The  flow  was  intermittent  between  0.4  <  y/6.  <  0.9 

w 

while  it  is  so  between  0.4  <y/<5j.  <1.2  in  the  natural  boundary  layer.  Thus 
in  the  synthetic  turbulent  boundary  layer  all  large  coherent  eddies  may  be 
considered  equal  while  they  vary  in  size  in  the  natural  case.  This  ine¬ 
quality  may  be  traced  to  the  various  origins  of  the  large  eddies  which  are 

probably  quite  randomly  located  in  the  turbulent  boundary  layer. 

The  engulftoent  of  non  turbulent  fluid  by  span wise  oriented  vortices 
was  considered.  Firstly  an  attempt  was  made  to  estimate  the  rate  of  entra¬ 
inment  on  a  plane  of  symmetry  of  an  isolated  spot  from  the 

ensemble-averaged  data.  Such  estimates  were  made  previously  by  CCD  and 

WSF;  the  novelty  in  the  present  approach  stems  from  the  following: 

(1)  Both  streamwis e  and  normal  velocity  components  in  the 

X-Y  plane  were  actually  measured,  while  in  previous  in¬ 
vestigations  the  normal  component  was  calculated  by  as¬ 
suming  two  dimensional  flow. 

(ii)  The  average  celerity  of  the  interface  was  determined  in 

the  X  and  Y  directions  from  a  large  number  of  measure¬ 
ments  at  successive  downstream  locations. 

Previously  the  celerity  of  the  interface  was  calculated  assuming  coni¬ 
cal  similarity  (CCD);  it  was  assumed  to  be  constant  across  the  boundary 
layer  (WSF).  The  calculations  based  on  the  recent  data  suggest  that  the 
leading  interface  of  the  isolated  spot  entrains  non  turbulent  fluid  much 
more  vigourously  than  previously  believed.  The  same  entrainment  oalcula- 
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tions  were  applied  to  the  array  of  interacting  spots  and  were  compared  to 
the  entrainment  required  to  sustain  the  overall  rate  of  growth  of  the 
spots.  It  transpires  that  the  two-dimensional  calculations  of  entrainment 
can  account  for  20$  of  the  total  rate  of  entrainment  necessary  to  maintain 
the  known  rate  of  growth,  the  rest  has  to  come  from  vortices  oriented  in 
the  streamwise  direction,  which  engulf  non  turbulent  fluid  in  the  X-Z 
plane.  Streamwise  vortices  are  clearly  observed  by  flow  visualization 
(CCD,  GBR,  Matsui,  1980  and  most  recently  Carlson  and  Widnall,  1981)  they 
are  also  observed  in  a  free  shear  layer  (Bridental,  1978)  and  are  largely 
responsible  for  the  mixing  of  species. 

The  use  of  a  hot  wire  rake  which  spans  across  the  boundary  layer  en¬ 
ables  one  to  consider  the  instantaneous  velocity  field  and  compare  it  to 
the  ensemble  average  data.  It  is  obvious  that  considerable  amount  of  in¬ 
formation  is  lost  by  the  averaging  process.  A  single  spot  contains  several 
large  coherent  eddies  and  an  array  of  successive  spots  after  a  prolonged 
interaction  contains  twice  as  many  eddies  as  spots  in  the  array.  The  iso¬ 
lated  spot  is  not  equivalent  to  a  single  eddy,  however,  the  eddies  within 
the  spot  must  be  arranged  in  some  preferred  order  to  give  the  spot  its  un¬ 
iversal  shape. 

Spots  generated  at  the  same  instant  at  different  spanwise  locations 
appear  to  grow  independently  of  one  another  and  not  inhibit  each  others 
growth  (see  also  Elder  I960).  Successive  spots  generated  from  the  same 
souroe  inhibit  the  growth  of  one  another  an  a  result  of  the  interactions. 
In  terms  of  ensemble  averaged  quantities,  this  inhibition  is  manifested  by 


Page  71 


an  acceleration  of  the  trailing  interfaces  in  an  array  of  spots  relative  to 
the  isolated  spot  surrotnded  entirely  by  laminar  flow.  One  ma.f  argue  that 
the  orientation  of  the  individual  eddies  in  a  spot  determines  the  effect  of 
its  interaction  with  surrounding  similar  structures.  Since  the  eddies  are 
basically  oriented  longitudinally  in  the  streamwlse  direction  the  interac¬ 
tion  of  successive  spots  may  severely  Inhibit  their  growth  process,  while 
interaction  with  neighbouring  spots  in  the  spanwise  direction  represents  no 
restriction  on  the  growth  of  these  eddies. 

The  conventional  statistical  methods  used  presently  to  describe  the 
mean  flow  in  a  transitional  spot  are  inadequate.  An  attempt  to  present 
physioally  meaningful  ensemble  averaged  turbulent  intensities  in  this  flow 
(see  Antonia,  Chambers,  Sokolov  and  Van  Atta  1981)  is  even  less  adequate, 
beoause  it  incorporates  errors  resulting  from  variations  in  the  time  of  ar¬ 
rival  of  the  spot  at  the  measuring  location;  jitter  in  the  length  of  the 
spot  etc.  Some  of  these  effeots  have  been  diagnosed  and  assessed  but  there 
is  no  dear  alternative  description.  One  should  define  a  statistical  re¬ 
presentation  which  will  reveal  the  detailed  structure  of  the  ensemble  aver¬ 
aged  spot  and  at  the  same  time  ensure  that  it  is  representative  of  the  sin¬ 
gle  realization. 
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FIND  THE  FOLLOWING  COEFFICIENTS: 

aq,  bo,  Ca.  Aq,  Bq  AND  Cq, 

OF  A  TWO  DIMENSIONAL  LINEAR 
POLINOMIAL  BY  SOLVING  SIX 
EQUATIONS  WITH  Q  AND  a  KNOWN 
AT  EACH  POINT. 

THREE  EQUATIONS  FOR  Q 

Q(I+1,J)  -  Q(I.J)  +  DE*Bp 

Q(I,J+1)  *  Q(I,J)  +  DE*Aq 

Q(  1+1  ,J+1)  *  Q(I,J)  +  DE*Ap  +  DE*Bp  ♦  DE2*C( 

THREE  EQUATIONS  FOR  a 

ot(I+1,J)  -  a(I,J)  +  DE*Ba 

a(I,0+1)  -  a(IfJ)  +  DE*Aa 

o(I+l.J+1)  -  o(I,J)  +  DE*Ajj  +  DE*Ba  +  DE2* 


PLAY  BACK  EACH  X-WIRE  DATA  AND 
FIND  THE  NARROWEST  WINDOW  WITHIN 
THE  CALIBRATION  REGION  CONTAINING 
ALL  DATA  POINTS. 
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STEP  7 


VOLTAGE  TO  VELOCITY  TRANSFORMATION 


FIND  THE  SQUARE  IN  THE  ABOVE  WINDOW, 
CONTAINING  EACH  MEASURING  POINT  e2- 

e1 

I  -  f(fif) 


J  -  f<n|> 


Pag-  75 


I  I4I,J  ♦  I 

- 1 - * 


T  XJ1 

DE2 


\hJ+\ 


H-DEI 


STEP  8 


TRANSLATE  VOLTAGES  TO  VELOCITItS 
AND  ANGLES,  USING  TWO  DIMENSIONAL 
LINEAR  POLINOMIAL 


Q  =  Vj+Vde1+bq^E2+Vde,*DE2 

a  =  «i  ,+A  *DEUB  *DEZ+C  *DE1  *PE2 

*  U  (I  (J 


STEP  9 


TRANSFORM  COORDINATES 
U  =  Q*COSot:  V  =  Q*SIN« 
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LIST  OF  SYMBOLS 


A,B,C,D, 


-  Calibration  constants  of  normal  hot-wires. 


A  ,B  ,C 

a  a  a 

VVCQ 


El  ,E2 


-  Calibration  constants  of  X-wires. 


-  Hot-wire  voltage. 

-  Voltages  of  X-wire. 

-  The  mean  height  of  a  spot  in  a  laminar  boundary 
layer . 


I  -  Intermitteney  function. 

o 

(2/pU  )x(dp/dx)  -  Pressure  gradient  along  the  flat  plate. 

Q  -  Resultant  velocity  during  X-wire  calibration, 

q  -  Rate  of  entrainment  (velocity  times  length). 

T  -  Time. 

AT  -  Time  interval  between  successive  sparks. 


TnsATxU  /X. 

oo  s 

U.V.H 


-  Dimensionless  interval  between  sparks. 

-  Components  of  instantaneous  velocity. 


u,v,w 

OaU-U 


-  Components  of  velocity  fluctuations. 

-  Mean  value  of  streamwise  velocity  component. 


Uint.Vint 


-  Streamwise  and  normal  components  of  the  celerity 


Oi,»i 


Vi. 

V’p 


of  the  interface. 


-  Mean  value  of  streamwise  and  normal  velocity  in  the 
laminar  neighbourhood. 

-  Two  components  of  the  celerity  of  the  leading  edge. 

-  Tw©  components  of  perturbation  velocity. 
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UT  -  Friction  velocity. 

U/U  =U+  -  Dimensionless  velocity. 

T 

<U>  -  Ensemble  averaged  velocity, 

u’ ,v'  -  R.M.S  of  u  and  v  fluctuations. 

<u'>,<v»>  -  R.e.m.s  (root  ensemble  mean  square)  of  the 

fluctuations. 

U’v'  -  Mean  product  of  uv. 

<uv>  -  Ensemble-averaged  of  uv  product. 


u*  ,v* 

7SS 

u*  v’ 

X.Y.Z 

x,y,z 

Xs 

yUT/v=y+ 

a 


-  Turbulent  zone-averaged  of  the  fluctuations  u  and  v. 

-  Turbulent  zone-averaged  of  uv  product. 

-  Free  stream  velocity. 

-  Directions  in  Cartesian  coordinates. 

-  Distances  along  three  directions. 

-  Streamwise  distance  from  the  spark  generator. 

-  Normal  distance  from  the  surface  in  wall  coordinates. 

-  Direction  of  flow  during  X-wire  calibration. 

-  Intermitteney . 

-  Boundary  layer  thickness. 

-  Laminar  boundary  layer  thickness. 

-  Turbulent  boundary  layer  thickness. 

-  Displacement  thickness. 

-  Momentum  thickness. 


Red*  -  Reynolds  number  based  on  displacement  thickness. 

Re6  -  Reynolds  number  based  on  momentum  thickness. 
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